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A MODEL FOR SIMULATING AN AIR-TRANSPORTATION SYSTEM* 


D. A. D'Esopo, H. L. Dixon, and B. Lefkowitz 


Stanford Research Institute 
Menlo Park, California 





INTRODUCTION 

This paper is concerned with a model for simulating the movement of cargo over a 
complex transportation network that may contain transloading points and in which several dif- 
ferent types of transport aircraft may be used. The primary quantities of interest are the 
number of transport aircraft of each type that are required in each phase of a phased airlift 
operation and the duration of the start-up phase of an airlift. The start-up phase is considered 
to be complete when the subsequent flow of cargo to the final destination can be accurately 
approximated by a rate of delivery. Since normally this is the case after items designated as 
high priority have arrived, the duration of Phase 1 can be thought of as a measure of the ability 
of the transportation system to complete the urgent part of a new assignment. The scheduling 
of transportation operations in such a network bears some analogy to industrial production 
scheduling, and it is possible that the method of approach described here may be applied to 
problems in the latter area [1-4]. 

It may be convenient here to interpret the word 'model" as used in the title of this 
paper, Underlying the hand calculation of the numerical consequences of an airlift—e.g., time 
to complete a mission —is some conception of how the consequences depend on both the deci- 
sions about the handling of the air-transport operation and such things as the performance 
characteristics of the aircraft. When the researcher is attempting to determine the conse- 
quences in any specific case he will make whatever reasonable assumptions are required to 
get a reasonable answer from whatever description of the problem is at hand. In other words, 
his assumptions and technique for solving the problem can be tailored to each case. When, 
however, this task is to be referred to an electronic computer, these assumptions have to be 
made beforehand, and they have to be reasonable for a sufficiently wide range of problems to 
justify the labor of coding the computations. Moreover, the format of the description of the 
problem must be set forth beforehand and, again, must be flexible enough so that a reasonably 
large number of problems can be described—and yet not be so cumbersome that every problem 
tas tobe worked dozens of times. It is with both the format for statement of the problem and 
the assumptions underlying the determinations of the results that the term "model" is 
concerned, 

While it is conceivable that a model of the start-up operation of an airlift might take 
the form of a set of equations relating the quantities of interest to the pertinent system vari- 
ables, it did not appear possible to formulate a model of this type which would be applicable to 
Miiitiet received July 13, 1959. 
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a sufficiently large set of cases. Consequently, a different type of model was developed: one 
by which the performance of the real system during this start-up period can be simulated, 

The development of such a model consists of two steps. First, a set of parameters 
must be determined that are sufficient to describe the state of the system at any point in time, 
These parameters are the time and place of availability and the physical characteristics of 
both cargo and transport aircraft. Second, a set of rules must be specified that determine the 
transformation of the parameters during each successive interval of time. The rules, which 
could be either stochastic or deterministic, dictate cargo assignments and transport-aircraft 
movements, Repeated application of these rules to the parameters will generate a sequence of 
states, This sequence of states is intended to approximate the actual sequence of states of the 
real transportation system. If now the rules are programmed for a computer and the initial 
state of the system is described, the machine can simulate the entire operation. From the 
simulation the number of transport aircraft, the duration of the operation, and a variety of 
other information can be obtained. 

In the organization of this paper the parameters which describe cargo, transport air- 
craft, and depots are examined first. (Depots, as used here, are all locations where aircraft 
pick up or deliver cargo.) Then the rules for transforming the parameters are discussed; that 
is, the scheduling rules. Next, the interaction between Phase 1 and Phase 2 of the airlift is 
examined. Then the method by which the model is applied to the problem of finding the number 
of carriers needed to complete Phase 1 in the shortest time is considered. Finally, the way 
transport aircraft requirements for the Phase 2 operation can be handled analytically is exam- 
ined briefly. 

The language used in this paper to describe the model is chosen so as to be readily 
translatable into computer operations. For example, the word "table" is used instead of the 
more abstract word "vector." 


PARAMETERS 
Cargo 
The amount of cargo that can be carried by one transport aircraft in a single shuttling 


movement! is limited by either the cubic capacity of the aircraft or the allowable cargo weight. 


In this model, cargo is described in terms of weight rather than volume. 

A homogeneous treatment of cargo weight, however, would lead to certain erroneous 
results, as, for example, the assignment of packaged cargo or personnel to tanker aircraft. 
For this reason, independent schedules are drawn up for each cargo class. Furthermore, to 
insure realistic schedules it is necessary to distinguish the high-priority cargo, called Phase 
1 cargo, from lower-priority cargo, called Phase 2 cargo. High-priority items are those that 
must be moved as soon as transport aircraft and cargo availability permit. 

In addition to the cargo class and priority, the initial geographical location of the items 
must be specified. With this information, it is possible for the model to develop what is called 
the Weight-at-Depot Table. This table shows the initial disposition of pounds of items to be 
moved, indexed by their location, cargo class, and priority. To facilitate the creation of this 
table, a code (the prepositioning class code) is used to identify items having similar storage, 
cost, weight, or essentiality characteristics. This code greatly simplifies the specification 
and revision of cargo prepositioning strategies. In particular, the initial location of items is 








ISee Transport Aircraft section. 
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specified by statements of the form, "such-and-such percent of so-and-so cargo class is to be 
located at depot X,"" and these statements are combined in an obvious fashion with detailed 
descriptions of line items to produce the aforementioned table. 


Transport Aircraft 

In the model, the logic for routing the transport aircraft assumes that the aircraft 
engage in shuttling movements, i.e., they always return to their original loading point for pos- 
sible reloading after they have been unloaded at their destination. By restricting aircraft 
routing strategies to this one mode of operation, the requirement for computer storage is 
greatly reduced. If nonshuttling movements are allowed, then it would be possible to dispatch 
any available aircraft to any depot in the network. To keep track of such a network would 
require computer storage of: (1) the status of all inbound shipments at every depot and (2) the 
performance of each aircraft type over every admissible route. In addition, if nonshuttling 
movements were allowed, it would be necessary to develop a procedure by which the computer 
selects destinations and routes whenever an aircraft is scheduled. Conceivably, this selection 
could be based on the outcome of a series of experiments conducted within the computer, or on 
the adoption of actual route-selection procedures, or on the application of some optimal rule. 
However, using a series of experiments would probably necessitate an excessive amount of 
computing time; the adoption of actual route-selection procedures presumes the existence of 
procedures that are flexible enough to cover all situations, a state which does not seem to 
exist; and the specification of an optimal rule escaped the ingenuity of the authors. 

By restricting transport aircraft movements to shuttling movements an important sim- 
plification occurs: the entire schedule of aircraft movements can be obtained one depot at a 
time, provided that depots are treated in an order compatible with the shuttling strategy. Thus, 
the situation is reduced to one in which it is necessary to know the status of inbound shipments 
at the particular loading point being analyzed (and at no other). Furthermore, for any aircraft 
in the system it is necessary to know its performance only over the route it travels. 

It should be noted here that shuttling movements of the type just described cannot handle 
cases where aircraft loads can be adjusted for different legs in a route. An example is illus- 
trated in the following sketch. 


Capacity = 15,000 pounds Capacity = 30,000 pounds 


— ie 
New York — 


pte sudh eo nhaaia geen San Francisco 


Capacity = 10,000 pounds 


Suppose a single aircraft is required to move 30,000 pounds of cargo from New York to 
San Francisco and only one intermediate stop, Denver, can be used. Suppose that in a single 
movement the aircraft can carry 10,000 pounds in a nonstop flight from New York to San Fran- 
cisco, 15,000 pounds from New York to Denver, and 30,000 pounds from Denver to San Fran- 
tisco, One possible method to move the 30,000 pounds is to make three nonstop trips from 
New York to San Francisco carrying 10,000 pounds each trip. However, if loading and unload- 
ing time represents only a small portion of total trip time, the length of time to move all the 
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material to San Francisco is reduced by moving 15,000 pounds from New York to Denver, 
returning to New York to move the second 15,000 pounds, stopping at Denver to pick up the 
15,000 moved previously, and proceeding to San Francisco with the 30,000 pounds, If only 
shuttling movements are allowed, then in the situation just described, the computer would have 
to consider the aircraft flying from New York to Denver to be different from the one flying 
from Denver to San Francisco. Despite inadequacies in situations of this type, shuttling move. 
ments will cover scheduling strategies of great interest to planners of air transportation, 

In the computation system, the inbound and outbound portions of the shuttle are permit. 

ted to include refueling stops, and the unloading point need not be the ultimate destination of 
the operation; that is, the system will handle cargo transloading. Further, aircraft do not have 
to be refueled at the final destination, but can fly a radius mission from their last fueling point, 

The two main consequences of restricting aircraft to shuttling movements are: 

1. Each aircraft can be represented by a small number of parameters, i.e., origin, 
destination, travel time, fuel consumption, and cargo capacity. The cargo capacity 
for a particular aircraft-route combination depends on the physical parameters of 
the aircraft, the length of the longest leg of the route, and a decision as to whether 
or not the aircraft is refueled at the destination. 

. Aircraft possessing the same parameters can be aggregated into groups, and certain 
group-wide parameters can then be introduced. (A group is the set of aircraft, all of 
the same type, that originates at the same airfield, travels the same route, and 
carries the same class of cargo.) 


The computer stores all the pertinent data concerning the characteristics of aircraft ina 
Group-Parameter Table. 





The computer has a third table—the Aircraft-Availability Table—which contains the 
status of all aircraft to be loaded at the particular depot under analysis. Each aircraft is 
tagged by its group identification and time of availability. Initially, the time of availability 
represents the time the aircraft can first be scheduled. By cross references to the Group- 


Parameter Table, it is possible to determine the performance characteristics of each aircraft 
as it becomes available for scheduling. 





Depots 

As stated earlier, the entire schedule of aircraft movement in the network can be 
obtained one depot at a time, provided that depots are treated in a particular order. The order 
is correct if no depot is treated before all its suppliers have been processed. Such orderings 
are not necessarily unique, but once the shuttling strategy is specified, the permissible order- 
ings are readily constructed. The order of treatment is stored in a table called the Depot- 
Order Table. 

At any point in the analysis, the status of all outstanding shipments in the network is 
contained in a table called the Unprocessed-Shipments Table. Each shipment is described by 
weight, destination, and time of arrival. The initial state of this table is obtained from the 
Weight-at-Depot Table by regarding the Phase 1 cargo initially located at a depot as a ship- 
ment that arrived at time zero. After establishing the identity of the next depot in the analysis, 
through reference to the Depot-Order Table, all inbound shipments to that depot can be removed 


from the Unprocessed-Shipments Table and stored in increasing time sequence in a final table, 
called the Inbound-Shipments Table. 
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SCHEDULING RULES 

At this point all the necessary preliminaries have been performed. That is, the com- 
puter knows how much cargo is to be moved, where it is located, and when it will be ready for 
movement, The computer also knows where the transport aircraft are located, what their per- 
formance characteristics are, and when they will become available. Finally, the computer 
knows the routes, the final destination, and the order in which the depots are analyzed, It is 
now possible for the computer to schedule shipments from a depot. Basically, this consists of 
a series of decisions and operations on successive entries in the Inbound-Shipments Table and 
Aircraft-Availability Table, to produce new entries in the Unprocessed-Shipments Table. 

The schedule of shipments from a depot is obtained by the repeated application of the 

following steps: 

1. Determine from the Aircraft-Availability Table the time of availability and group 
identification for the next available aircraft. 

2. Determine from the Inbound-Shipments Table the time at which the next aircraft 
load becomes available.” 

3. Record the shipment size, destination depot, and time of arrival in the Unprocessed- 
Shipments Table. The time of arrival equals departure time plus outbound travel 
time plus adjustments for loading, unloading, and en-route stops. 

4, Finally, modify the Available-Aircraft Table to show the time when the aircraft just 
scheduled will be next available. This equals the departure time plus the round-trip 
travel time, including the loading, unloading, and en-route stop times. 

The scheduling of aircraft from a depot ends when the total weight of scheduled shipments 
equals the total weight of all entries in the Inbound-Shipments Table. That is, there is nothing 
more to be moved. The above procedure is applied to each depot in the network in the order 
specified by the Depot-Order Table. 


SUPPLIES IN PHASE 1 

Thus far the analysis of Phase 1 has assumed that the total weight to be moved is com- 
pletely specified in the initial setting of the Weight-at-Depot Table. However, some of the 
wits moved during Phase 1 might start consuming supplies immediately upon arrival at their 
final destination. Therefore, enough supplies (Phase 2 cargo) must be carried in Phase 1 to 
support these units until the first deliveries of resupply arrive in Phase 2. Three important 
elements of data must be determined before the amount of these initial supplies can be com- 
puted, They are: (1) the schedule of arrivals of supply-consuming units at the final destination, 
(2) the time that resupply begins to arrive at the final destination (i.e., the time when the 
earliest of the Phase 2 shipments arrives at the final destination), and (3) the rates at which 
supplies are consumed. The latter is specified in the input data. 

It should be noted here that the amount of time needed to complete the Phase 1 opera- 
tion cannot be determined until the amount of supplies that must be moved in Phase 1 is known. 
But the amount of supplies that must be moved in Phase 1 depends on the amount of time needed 
for the completion of Phase 1, which in turn depends on the amount of Phase 1 supplies. To 
break this apparent impasse, a scheme of successive approximations is used. First, complete 
schedules are drawn up for the deployment of Phase 1 cargo with no supplies. Then estimates 
a 


d , 
The aircraft departs at the later of these two times. The load it carries is either (a) its 


‘apacity over the route it must travel or (b) the total of the remaining cargo to be shipped, 
chever is smaller. 
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of supply consumption are made on the basis of this tentative schedule. These estimated Phay 
1 supplies are then apportioned to the various depots in accordance with a prepositioning stray. 
egy, and the scheduling process is repeated with the new weights. If desired, the amount of 
Phase 1 supplies can be re-estimated on the basis of this new schedule and the scheduling 
repeated a third time. 


MINIMUM NUMBER OF CARRIERS IN SHORTEST TIME 

The program user can specify the number of transport aircraft that are expected to be 
available, and the model will use only those aircraft in developing schedules, Alternatively, 
the model can compute the minimum number of aircraft needed to complete Phase 1 in the 
shortest time. In either case, the entire simulation can be performed for different numbers of 
aircraft in each group. 

If group sizes are varied, it seems obvious that beyond a certain point increasing grow 
sizes will not decrease the completion time for the entire operation; that is, the operation will 
require at least a one-way flight time, regardless of fleet size. A series of experiments could 
be performed to find the smallest group sizes that complete the operation in the shortest time, 
The results of this experimentation are especially meaningful when the distribution, or mix, of 
the aircraft types at each loading point is fixed. If the mix at an airfield were not specified, 
then the computer would always choose the most productive aircraft available to it, and the 
task of finding the minimum numbers of these aircraft would be trivial. The computer does 
perform automatically a series of experiments to find the minimum numbers for each group, 
In particular, the computer will find the minimum numbers of aircraft required to obtain either 
one of two "shortest time" objectives: 

1. To complete Phase 1 movement from each depot in the shortest time permitted by 

the scheduling procedure and aircraft capabilities, or 

- To complete Phase 1 movement from each depot in the greater of the time just men- 
tioned or a time that will permit the last aircraft departing from a depot to reach its 
destination no later than the time that any aircraft operating from any previously 
analyzed depot reaches its destination. 

In order to use either of these options, the composition (mix) of the aircraft fleet at 
each airfield must be specified. That is, the distribution among groups of the total number of 
aircraft required to move a given cargo class from a given airfield must be known. The dis- 
tribution is specified by stating the percent of the total required aircraft at a depot that belong 
to each group assigned there. This percentage is called the fleet mix percentage. 

Since the fleet mix percentage must be specified, it is more precise to state that the 
model will determine the number of aircraft in each group required to obtain either of the 
above two shortest time objectives in such a way that the following conditions are met: (1) the 
aircraft at any airfield are distributed among the groups in proportion to the fleet mix percetl- 
ages and (2) the total number of aircraft at each airfield is as small as possible for the fleet 
mix percentages specified. 

When the first time objective is chosen (i.e., when it is desired to complete Phase 1 in 
the shortest possible time) the critical time used in the trial-scheduling procedure is taken to 
be either the time the last shipment of cargo arrives at the airfield or the time the last 
required aircraft is available, whichever is later. 

Objective 2 attempts to "balance" time requirements throughout the network. For 
example, suppose it is desired to airlift equal amounts of cargo from San Francisco and 
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Washington, D. C., to New York, and that all cargo from both San Francisco and Washington 
must be delivered to New York before the airlift can be considered complete. Suppose, too, 
that all of the required aircraft are immediately available. If in this situation Time Objective 
i was specified, the resultant fleet at Washington would be larger than necessary. Figure 1 
ilustrates this point. 

In Figure 1 the curves show the relationship 
between time to complete Phase 1 and the number of air- 
craft needed for both the San Francisco (SF) and Washing- 
ton (DC) operation. If Time Objective 1 is specified, the 
umber of aircraft required and the time to complete the 
iirlift are indicated by Nj, T; for Washington and No, Tg 
for San Francisco. If Time Objective 2 is specified, the 
Washington-based fleet must be capable of completing its 
part of the airlift no later than the landing time of the last 
San Francisco-based aircraft. By referring again to 
Figure 1, it can be seen that, given Ty to complete its job, 
the Washington-based fleet can be as smallas Ng aircraft. 
Both time objectives result in the same number of aircraft 
and the same completion time for the San Francisco lift. 
Thus, if Time Objective 2 were specified instead of Time pen SF see 
Objective 1, a total of (N, - Ng) fewer aircraft would be Figure 1 
needed, even though the entire lift would still be completed 
by time Ty. 

As before, the computer searches for the minimum aircraft fleet size at each depot, one 
depot at a time, in the order specified by the Depot-Order Table. The technique used to obtain 
the minimum total aircraft fleet size at any depot is as follows. First, determine the comple- 
tion time for that depot according to the stated time objective. When Time Objective 1 is cho- 
sen (i.e,, when it is desired to complete Phase 1 in the shortest possible time), then the com- 
pletion time for the depot equals the time required (e.g., T, and Ty in Figure 1) by an aircraft 
fleet of the correct mix whose aggregate capacity equals the total Phase 1 weight that must be 
lifted from that depot. This completion time is obtained by actually scheduling individual air- 
craft movements for just such a fleet. When Time Objective 2 is specified, the completion 
time at the depot is the later of the times just described and a time equal to: 


TIME TO COMPLETE PHASE I 











The maximum of the arrival times of all shipments scheduled from previous depots 
minus 
The maximum of the times needed by groups to deliver cargo from the depot.> 


After the computer calculates the lift-completion time for a depot, it prepares a 
sequence of schedules for a succession of fleet sizes, increasing or decreasing the fleet size 
sed in the next trial schedule, depending on whether the previous schedule was completed 
ter or before the lift-completion time. 


a 


iy, ‘ : — , , 
Note that this statement applies to the origination depot, but the completion times used for 
illustration in Figure 1 apply to the destination depot. 
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PHASES 2, 3, AND 4 CARRIER REQUIREMENTS 

Up to three additional airlift phases can be handled by the system. The basic differ. 
ences between Phase 1 and these last three airlift phases arise from the fact that the duration 
of the last three phases is specified and that cargo is assumed to be delivered at a uniform 
rate (e.g., so many pounds per day) for the duration of each phase. Recall that, in the simula- 
tion, Phase 1 cargo was scheduled through the system in individual aircraft loads. Since 
delivery rates are being considered in the last three phases, the computation method for deter. 
mining the number of aircraft is much more direct and simple. One additional assumption is 
needed to make the computations valid: the same amount of time is available at each of the 
depots in the network for the completion of any one of the last three phases. 

The number of transport aircraft required to complete each of the last three phases 
can readily be computed by using the transport aircraft utilization rate, the amount of cargo to 
be delivered, and the data in Phase 1 describing transport-aircraft performance over the 
selected routes. 

Phase 1 and the subsequent phases can be used together or independently in simulating 
an air-deployment operation. For example, if urgency of completion is of paramount impor- 
tance, all items to be moved could be moved in a Phase 1 type of operation. If a reasonably 
steady state of delivery is desired, only Phases 2, 3, and 4 need be considered for the analysis, 
If a portion of the items to be moved must be delivered as soon as possible and others can be 


delivered at some later time, then Phase 1 and the subsequent phases of the system would be 
used. 
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SUMMARY OF A METHOD FOR 
DETERMINING THE MILITARY WORTH OF SPARE PARTS*t 


Marvin Denicoff, Joseph P. Fennell, and Henry Solomon 


The George Washington University 
Logistics Research Project 


INTRODUCTION 


The principal purpose of this paper is to describe a method which has been developed 
for determining the military worth of spare parts, an approach which is currently being imple- 
mented in the Navy. The general frame of reference is the need of military worth information 
for inventory problems. This necessity is obvious since the main reasons for holding inven- 
tories are to support some current strategic plan and to present the feasibility of future and 
ucertain strategic plans. The specific frame of reference is the allowance list or shipboard- 
inventory problem. 

While it has been recognized for some time that some concept and measurement of 
military worth is necessary for inventory systems, this need became pressing in the light of 
observed characteristics of demands for spare parts. The results obtained from an analysis 
of usage data for 12 submarines over a four year period [1] pointed to the immediate necessity 
for developing a feasible and operational approach to military worth. 

The requirements study revealed that approximately 75 percent of the installed techni- 
cal items (the population range) showed zero movement by or on account of each individual sub- 
marine over a four-year period. Of the items which did move, more than 70 percent were 
demanded only one time over these four years. The problem was further complicated by the 
lack of repetitive movement; for any one ship, less than 1 percent of the technical items used 
inthe four-year period were used in every.one of the four years. Since this study, other 
empirical investigations have yielded the same results. 

The problem, at least a major aspect of the allowance list problem, thus resolves itself 
into the development of a method for treating those items which have shown zero usage in the 
past but which may move in the future. More particularly, we are concerned with pinpointing 
that segment of the population range which is vital to the ship's existence, independent of the 
mst history of usage. Shipboard stowage constraints (most critical in the case of submarines) 


und budgeting considerations are factors which influence this concern in the face of the usage 
fattern just described. 


SD 

‘Manuscript received April 23, 1958. 
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This paper describes a technique for making military worth evaluations. It discusses 
questionnaires developed for this purpose and provides details on a program which represent; 
a first attempt at achieving what is basically a relative ranking of the importance of repair 
parts equated to the essentiality of parts to components and of components to the mission of 
the vessel. Some 31,600 repair part applications (parts tied to parent components) are includ 
in the study. 

In an earlier paper in this Journal [2], Karr has described an approach which is similx 
in certain respects to that being described here. However, there also are some important dif- 
ferences in the definition of "worth," the procedures followed, and the use of the results. 


THE APPROACH 

One of the goals of our research, predicated on CNO criteria, is to develop Allowance 
List techniques which maximize the independent afloat endurance of combatant vessels. Inde- 
pendent operation, in this case, refers to the capability of ships for accomplishing tactical 
missions independent of reliance on shore base, tender, repair ship, or fleet-issue ship supply 
support. The Navy's retaliatory power in the event of enemy attack is, to some extent, depeni- 
ent on the degree of self-sufficiency which is achieved for its fighting ships. 

In this context, military worth may be thought of as a relative ranking system which 
measures the effects of shortages on the ship's capabilities. There are, in a general sense, 
two factors which influence the seriousness of a particular shortage. One of these is mission 
effect, which measures the effects of the inoperability of specific components on the ship's 
capability for accomplishing its assigned mission. The other, maintenance potential, has to do 
with the effect of end item or part failures on the operability of the parent component. Where 
such failures would render the parent component inoperable, maintenance potential considers 
the capability of the ship's force, in the event of a part shortage, to maintain the component in 
a satisfactory operable condition through on-board manufacture of the required part, and/or 
cannibalization, and/or the employment of jury-rigging procedures. 

The combining of these two elements, mission effect and maintenance potential, permits 
a relative evaluation of the seriousness of repair part shortages. For example, we can con- 
ceive of an item worth or shortage penalty scale which includes extremes as represented in 
Chart I. 

Determinations of spare part essentiality, it can be seen, are a combination of both of 
these elements — mission effect and maintenance potential. This is most clearly evident when 
we think in terms of parts with similar design characteristics or even the same part with 
multiple equipment applications. The worth of two similar parts will vary considerably if one 
is installed in a high-mission-effect component, while the other is related to a component with 
a low-mission effect. In the same manner, the worth of similar parts related to different com- 
ponents having equal mission effects would vary if the maintenance potential were high for one 
of these applications and low for the other application. 

There are, in fact, two separate decision processes involved in making spare parts 
military worth estimates. One, a determination of mission effect, is logically a command 
decision. The other, a determination of maintenance potential, is the decision of qualified 
maintenance personnel with the appropriate technical skills. 

This distinction had an effect on the preparation of the military worth questionnaires, 
and on the solicitation of personnel to participate in the program. 
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DETERMINING THE MILITARY WORTH OF SPARE PARTS 


CHART I 
Two Evaluations of the Seriousness of Repair-Part Shortages 





[ _ High Military Worth 


Low Military Worth 








——— 


Part a for Component W 


Mission Effect 





Inoperability of the component would 
necessitate termination of the mission. 


Maintenance Potential 





Part failure renders parent component 
inoperable. Required part not available. 
On-board manufacture, cannibalization, and 
jury-rigging procedures infeasible. 


Part b for Component X 


Mission Effect 





Inoperability of the component would have 
a negligible effect on the accomplishment of 
the mission. 


Maintenance Potential 





The parent component will continue to 
operate in a satisfactory condition without the 
necessity of replacing the failed part. Should 
replacement be desired, on-board manufac- 


ture, and/or cannibalization, and/or jury- 
rigging procedures are feasible. 














THE QUESTIONNAIRES 
Background 

Military worth questionnaires tailored to the separate elements of mission effect and 
maintenance potential were developed in August 1957. Since the immediate concern was the 
use of worth evaluations in the allowance list area, a decision was made to obtain worth esti- 
mates for the total component and part population range of a single combatant ship. The sub- 
marine, USS TIRU (SS 416), was chosen for this purpose. Population information for TIRU, the 


range and quantities of installed components and parts, was obtained from three Supply Demand 
Control Points. 


Mission Effect Questionnaire 

Description 

The questionnaire (Chart II) includes a description of a typical war-time patrol for a 
fleet-type submarine of the TIRU class. The patrol situation viewed by the questionnaire 
jarticipant as a set of assumptions which effect his evaluations, is as follows: 

(1) The submarine is to go on a patrol in a specified area for a period of 60 days. The 
mission of this patrol is to seek out any enemy shipping and sink same. In addition, the sub- 
harine is to report any incidental intelligence. 

(2) While on patrol duty, because of the area of operations, it is expected that the: sub- 
marine will be submerged an average of 18 hours per day. The submarine will be snorkeling 
he necessary amount of time. 

(3) It is to be expected that, during the 60-day period of operations, no supply replenish- 
nent can be made for installed items. In addition, no repair support will be available either 
{rom tenders, repair ships, or shore activities. 
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CHART 1 
Mission-Effect Questionnaire 





Mission Statement 





. Fleet-type submarine, TIRU class. 

. Sixty-day wartime patrol. 

. Submerged 18 hours a day. Normal snorkeling. 

. Complete isolation from supply or maintenance support. 





Assumptions 


-. In performance of the stated mission, the component under consid- 
eration fails. 
. Repair cannot be accomplished. 


Questions 


. Inoperability of the component would necessitate the termination of | 
the patrol action. 


- Inoperability of the component would introduce a high-risk factor in | 
the accomplishment of the mission. | 


. Inoperability of the component would introduce a moderate-risk factor | 
in the accomplishment of the mission. 


F 
O . Inoperability of the component would have a negligible effect on the | 
accomplishment of the mission. 





‘i 





Answering the Mission-Effect (Component) Questionnaire 

Agreement was reached to obtain three independent evaluations for each component. 
The questionnaires were divided among submarine officers with specializations in mechanical 
and electrical, electronic, or ordnance equipments. Nine men participated, three (for the 
three independent evaluations) assigned to each of the equipment categories. All of the par- 
ticipants were qualified submarine officers who have been executive officers or commanding 
officers on submarines. 

Component listings for TIRU were prepared from the population data decks supplied 
by the cognizant Supply Demand Control Points. These listings, one for each man, were pro- 
vided to the participating submarine officers as the media for recording the military-worth 
evaluations. 

The man answering the questionnaire assumes effectuation of the described patrol. He 
views the questionnaire from the perspective of the commanding officer of the submarine. 

The questionnaire describes four situations which may result from the failure of a 
particular component. These range from component failure necessitating termination of the 
patrol action, to failures which introduce high- or moderate-risk factors in the accomplishment 
of the mission, to failures which have a negligible effect. These situations or effects are 








into | 
capal 
as cl 
appl} 


rest! 
be cc 


cable 
4 pur 
poten 


avails 
These 


DETERMINING THE MILITARY WORTH OF SPARE PARTS 225 


designated by numeric evaluation codes: 1, 2, 3, or 4, respectively. The participant deter- 
mines which situation applies for each component in the study. Decisions are made by record- 
ing applicable code numbers to the component listings. 

Acode 1 decision (termination of patrol action) indicates a component failure, the 
seriousness of which would cause the ship to break off the patrol and immediately return to 
port for repairs. 

Acode 2 decision (high risk) indicates a failure which would introduce a calculated risk 
into the accomplishment of the mission, the risk being restrictive in terms of the operational 
capability of the ship. Depending on the type of component which has failed, limitations such 
as choice of areas of operation, selection of targets, reduced defense, capability, etc., might 
apply. The ship, however, would stay on station. 

A code 3 decision (moderate risk) indicates a failure which imposes a much less serious 
restriction on the accomplishment of the mission and wherein the component failure can often 
be compensated for; for example, by substitution of manual for mechanical operation. 

A code 4 decision (negligible effect) indicates a failure which imposes no restrictions 
on the accomplishment of the mission. 

In some cases where there are multiple quantities of particular components installed 
on the ship, the military worth evaluation varies with the quantity installed, e.g., there are 
four motors X installed on TIRU. If all four fail, the military worth evaluation is code 1; if 
three of the four fail, the military worth evaluation is code 2; if one or two fail, the military 
worth evaluation is code 4. Wherever such divisions were determined to be meaningful, the 
participant divided the total installed quantity into "significant splits" and provided a different 
aswer (code number) for each of these split quantities. This splitting factor, in a general 
sense, takes into account the "built-in" protection that is provided to the ship through stand-by 
equipments. 


Maintenance- Potential Questionnaire 

Description 

The questionnaire includes the same description of the wartime patrol that was appli- 
table to the mission-effect questionnaire. Again, the participants' evaluations are affected in 
apurposeful fashion by the assumptions inherent in the patrol situation. The maintenance- 
potential questionnaire is described in Chart II. 


Answering the Maintenance- Potential (Parts) Questionnaire 

Spare-part listings for TIRU were prepared from the population data decks supplied by 
te SDCPs. These listings were furnished to participating experienced military and civilian 
personnel at the SDCPs, to serve as the media for recording answers to the maintenance- 
potential questionnaire. Listings were prepared in stock number within component sequence 


so that parts with multiple applications would be listed as many times as there were applica- 
tons for the part. 





Participants were provided with data on the facilities for manufacture and hand tools 
Wailable on the TIRU. Information on the availability of bulk materials also was supplied. 
These data are essential to decisions on repair, manufacture, and jury-rigging capability. 

The participant, in this case, views the questionnaire from the perspective of a tech- 
lician serving on board the ship. He makes his decisions in consideration and knowledge of 
te maintenance capabilities that are feasible within the confines of the ship itself, with 
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CHART I 
Maintenance Potential Questionnaire 





. Replacement can be accomplished by the ship's force. 


. The required replacement part can be manufactured by the ship's 
force with the machinery available on the ship. 


. The equipment can be made to function through known Jury-rigging | 
procedures. No substitute spare parts are required. 


. The equipment will continue to operate satisfactorily for the stipulated | 
60-day period without a part replacement being made and without the 
necessity for jury-rigging. If there is a reduced efficiency, such 
reduction is acceptable for the operation of the equipment. 


T t 
O O 


. The part can be removed intact (without damage to the part) from the | 
equipment under consideration. 








deference being given to limitations imposed by the machinery, tools, bulk materials, and per- 
sonnel skills immediately available. 

All five questions are answered for each part for each of its applications; i.e., a part 
related to two different components requires two sets of answers, one for each of its compo- 
nent applications. Questions are answered as true or false. 

The participant examines each line item in relationship to its component application 
since decisions on installation capability, jury-rigging possibilities, and cannibalization poter- 
tial may vary considerably for similar parts or even the same part, depending on its component 
tie. For example, part "x'' might be installable by the ship's force for component "A," whereas 
the inaccessibility of component ''B" makes the installation of the same part "x" infeasible by 
the ship's crew. 

These refinements require the participation of technicians with a functional as well as 
design knowledge of the spare parts under consideration. Technicians utilized instruction books 
and BuShips' and manufacturer's plans in answering the questionnaire. 

Question 1 has to do with the capability of the ship's force for making the installation of 
the required part. Parts within the ship's installation capability are candidates for loading. 
Parts outside this capability would not generally be loaded since an emergency requirement for 
replacing such a part could not be satisfied by the ship's force. 

The remaining questions, predicated on an assumption of a part failure with no spare 
part in stock, have to do with the possibility of employing alternative resources for keeping 
the parent component operative for the stipulated 60-day patrol period. 

Question 2 lists on-board manufacture of the required part as an alternative solution t0 
the shortage problem. Repair capability of a damaged part is implied in this question, but is 
not listed explicitly. The reason for this lies in the subjective elements that enter into a 
repair decision; i.e., the extent and nature of the damage must be known and described before 
such decisions can be rendered intelligently. The large variety of damage possibilities 
precluded the use of this question in an explicit manner. 
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Question 3 deals with the possibility of jury-rigging. Preliminary orientation sessions 
with participants stressed the fact that only known or standard jury-rigging techniques were to 
be considered in answering this question. This stress on the employment of standard tech- 
niques reflects an interest in excluding from the questionnaire such subjective values as a 
requirement for “extraordinary talent." 

Question 4 describes a situation wherein the component will continue to operate satis- 
factorily despite the failure of the spare part under consideration. This situation, because it 
does not impose the burdens of on-board manufacture, jury-rigging, or cannibalization on the 
ship's force, represents the least serious type of failure. 

Question 5 treats the possibilities for cannibalization as an alternative solution to the 





shortage problem. It supplies, however, only a partial answer to the two-pronged aspect of the 


cannibalization decision; i.e., it asks only if the part can be removed intact (without damage to 
the part) from its parent components. Since it does not consider whether or not such removal 
will immobilize a component whose military worth is comparable to the component requiring 
repair, this question cannot fully answer the cannibalization problem. The total answer is 
derived from a consideration of Question 5 related to the military worth evaluations assigned 
tocomponents via the mission-effect questionnaire. For an example, let us conceive of a 

part "a'’ with common application to four components, W, X, Y, and Z. A negative answer to 
Question 5 (the part cannot be removed intact) has eliminated component W as a candidate for 
camnibalization consideration. Affirmative answers to Question 5 for the application of part "a" 
tocomponents X, Y, and Z have established these as possible candidates. However, suppose 
components Y and Z are eliminated on the basis of their mutually high military worth evalua- 
tions. Then, component X, with a low-worth rating, survives as the sole candidate for a can- 
nibalization requirement of part "a." 

An affirmative answer to any one of the foregoing questions has the effect of mitigating 
the seriousness of a shortage. We could additionally think in terms of a hierarchy of shortage 
values wherein a given shortage is rated as more or less serious, depending upon the number 
df alternatives there are to solving the problem occasioned by the shortage. 


DESCRIPTION OF RESULTS 
The computations described in this section of the report were performed for three 
purposes: 
- To check the consistency and accuracy of the answers relative to the military-worth 
questionnaires. 
- To establish categories of military worth. 
- To determine the frequency distribution of spare parts within these categories. 


Valuation of Components ! 

Earlier in this paper we described each of the four possible answers for any one com- 
ponent. In this section we shall present the quantities of components which fall into each of the 
four "worth" categories. In addition, since the answers by the participants include a personal 
hdgment, the consistency of answers given for each component among the three individuals is 
interest and will be described. 


‘When referring to a ''component" in this paper, we are, for mechanical and electrical compo- 
tents, using the term to represent "component service.'' Each component is evaluated in 
telation to the service it performs. For example, if there are two identical components 
stalled and each performs a different service, there may be a different valuation for each 
‘omponent unit. 
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It should be borne in mind that the tabulations in the following discussions are based op 
a conservative estimate of component worth. Where more than one unit of a component is 
installed and each unit performs the same function, the participants provided "split" evalua- 
tions. For example, if there are 10 units of a component installed, the answer may have been 
as follows: If 8-10 units became inoperative, the worth would be "1"; if 5-7 units became inop- 
erative, the worth would be "2"; if 3-4 units became inoperative, the worth would be "3"; and if 
only 1-2 units became inoperative, the worth would be ''4." In cases such as this, the value of 
the component is based on its highest worth. For example, in the above case, the worth assign. 
ment to each unit of the component would be "1." 

There are three general categories of components. These are (1) mechanical and elec- 
trical, (2) electronics, and (3) ordnance. For each one of these categories, three independent 
answers were provided for each component. The summary results will be presented in terms 
of these general categories. For the sake of brevity, detailed results will be presented only 
for "mechanical and electrical'' components. 


Mechanical and Electrical Components 

There were 1190 mechanical and electrical components evaluated. These are divided 
into Types 1 and 2. Type 1 refers to installed components, and Type 2 refers to portable com- 
ponents and equipage (e.g., messing gear, tools, office equipments, etc.). Type 1 consisted of 
692 components; the components in Type 2 numbered 498. 


The number of components assigned to each "worth" category by each of the three 
participants is shown in Table 1. 





TABLE 1 
Number of Mechanical and Electrical Components Assigned to 
Each Worth Category by Individual Participants 





Worth Category 
Components 





4 








Type 1 43 101 


Type 2 0 0 
Total 43 101 





Type 1 30 | 187 
Type 2 6 33 
Total 36 | 220 





Type 1 46 | 172 
T. Type 2 2 10 
Total 48 | 182 





























*Number of components not evaluated. 


One important observation from the data presented in Table 1 is that relatively few 
components have been assigned a high worth by any one of the participants. For example, Ty, 
assigned the highest worth category to only 43 components out of a total of 1190. 
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Let us now look at the degree of consistency among the three answers for each com- 
ponent. This is shown in Table 2. 


TABLE 2 
Number of Mechanical and Electrical Components in Terms of 
Range of Answers Provided by Participants 
(Consistency) for Each Component 





Range (R) 
Component 





Total* 








Type 1 (installed) 189 367 686 
Type 2 (equipage) 156 302 35 496 
Total 345 669 160 1182 


























*Excludes components for which at least one answer was not 
given. 


In Table 2 the range refers to the quantitative difference between the highest and lowest 
score for each component. Thus, for a component where R= 0, the three answers were in 
complete agreement. Where R = 1, two of the three answers were in agreement and the third 
differed by one worth category (e.g., 1,1,2 or 3,4,4, etc.). The type of information appearing in 
Table 2 is of great interest, as it sheds light on the validity of the approach. 

For approximately 80 percent of the Type 1 components, there was at most one dis- 
agreement among the three answers, and this disagreement was only in terms of one worth 
category. In the case of Type 2 components this percentage is over 90 percent. 

It is also of interest to add that in the great majority of cases where R= 1 the dis- 
agreement was between worth categories 3 and 4. For example, in the case of Type 2 compo- 
nents, 292 of the 302 components with R= 1, all three answers consisted only of 3's and 4's. 

Of particular interest is the number of cases where R = 3, for here at least one answer 
was 1 and one other answer was 4. The number of such occurrences was extremely few. In 
the few cases where this did occur, the reason was apparent. For example, one component for 
vhich at least two answers were at both extremes was a fire extinguisher. In this case two 
uswers were 4 and one answer was 1. The participant who answered 1 evidently assumed the 
ative outbreak of a fire, while the other two participants viewed this component from the 
standpoint of the probability of fire in the area where the extinguisher was located. There are 
afew such safety components where the interpretations of the questions given to the partici- 
pants could vary, and thus leads to a case* where R = 3. 

The number of components for which R = 2 is of some significance in that it represents 
inmost cases a clear difference of opinion among the three participants. In about half the 
cases where R = 2, each of the three answers was different (e.g., 1,2,3), while in the remainder 
cases two of the answers were in agreement (e.g., 1,3,3). Also, for those components where 
R= 2, very rarely was any one of the answers 1. The number of cases where some significant 


a 
Most of these components do not have spare parts. 


2500 ~ 60 - 3 
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differences of opinion arose, considering the varied individual experiences of the participants, 
seems small enough as not to disallow the validity and acceptability of the approach. 


Consolidation of Answers 





In order to arrive at a single worth category for each component, it was necessary that 
the three independent answers be consolidated. There are several possible rules for accom- 
plishing this consolidation. One such rule is to simply use the arithmetic mean of the answers, 
In view of the observed nature of consistency among the three answers, the arithmetic mean is 
an acceptable measure. This is so because in most cases very little variance is present. The 
number of all components in each worth category arrived at by computing the arithmetic mean 
of the three answers and rounding off to the nearest integer is shown in Table 3. 


TABLE 3 
Number of All Components in Each Worth Category 
after Consolidation of Three Independent Answers 





Components 





Mechanical and | Mechanical and 
Electrical Electrical Electronics 
(Type 1) (Type 2) 








31 0 
7 























Valuation of Parts 

The procedure adopted for determining the worth of parts relative to components in 
which they are installed was described earlier. Five questions were answered for each part 
application. 

For purposes of this study, all items were evaluated in regard to component applica- 
tions, independent of the particular cognizance Inventory Manager. Thus, for example, the 
participants at the Ordnance Supply Office evaluated electronics and general stores parts 
installed on ordnance components as- well as ordnance items installed on ordnance components. 
In the same way, the participants at the Electronic Supply Office and Submarine Supply Office 
made decisions on parts under the supply control of other inventory managers. 

This section describes results for technical spare parts installed only in "mechanical 
and electrical" components. The distribution of these parts in terms of worth categories pre 
scribed by Plan 1 (see Chart IV) is shown in Table 4. 

The results here are similar to those found for parts installed in electronics and 
ordnance components; i.e., a high percentage of cases wherein the ship's force could not com- 
pensate for a part shortage through its own devices. It should be remembered, however, that 
we are dealing with a submarine, a vessel which has extremely limited facilities for repair 
work. In this light, it is interesting to note that, even with such restricted facilities, 12 per 
cent of the technical items can be manufactured and/or jury-rigged. The 12.8 percent of the 
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TABLE 4* 
Distribution of Parts for Mechanical and 
Electrical Components Assigned to Worth 
Categories Prescribed by Plan 1 





No. of Part 


Category Applications Percent 








7414 75.2 
1180 12.0 
1258 12.8 


9852 100.0 




















*There were 3902 Category D parts. These 
were not included in the table. 


CHART IV 
Maintenance- Potential Categories 





Plan 1 Plan 2 





: No possibility of compensation. : No possibility of compensation. 
: Compensation through on-board man- : Compensation through jury-rigging 
ufacture and/or jury-rigging. only. 
No need for compensation. : Compensation through on-board man- 
Part not installable by the ship's ufacture only. 
force. : Compensation through cannibalization 
only. 
Compensation through any two of the 
three compensating factors. 
Compensation through all three of the 
compensating factors. 
No need for compensation. 
Part not installable by the ship's 
force. 
A>B>C A>B>C>ODOEDF>GOH 














items which could fail without disturbing the operation of the parent component also is a fact 
vorth noting. Certainly, with space as a constraint, these items would have low priority for 
on-board loading considerations. The distribution of mechanical and electrical parts in terms 
of worth categories prescribed by Plan No. 2 (see Chart IV) is shown in Table 5. 

Tables 4 and 5 are somewhat interesting in that they focus attention on the effects that a 
‘onsideration of cannibalization has on the possibilities of compensating the part shortages. 
The cannibalization possibilities here have a remarkable effect on the Category A classifica- 


tion; i.e., a reduction in the number of cases where there is no possibility of compensating for 
apart shortage. 
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TABLE 5* Combining the Part and Component 
Distribution of Parts for Mechanical and Evaluations 


Electrical Components Assigned to Worth Earlier in this paper a concept of mij. 
Categories Prescribed by Plan 2 tary worth was discussed wherein the seri- 
ousness of repair-part shortages was equated 





Category Ho. of Bust Percent to two factors, mission effect and maintenance 


Applications potential. Mission effect, obtained from the 
component questionnaires, relates to the effec 
of component failure on the ship's capability 
for executing its assigned mission. Mainte- 
nance potential, obtained from the parts ques- 
tionnaires, considers the capability of the 
ship's force, in the event of a part shortage, to 
maintain the equipment in a satisfactory oper- 
able condition through on-board manufacture 
9852 100.0 of the required part, and/or cannibalization, 
and/or the employment of jury-rigging 
procedures. . 








4069 41.2 
477 4.8 
103 1.0 
214 2.2 

3714 37.8 

17 0.2 

1258 12.8 





viawyFAoaw Pp 

















*Parts which cannot be installed by ship's 
force (Category H) are excluded. 
In order to obtain relative measures of 


the seriousness of part shortages, it becomes 

necessary to combine the mission-effect (component) and maintenance-potential (part) worth 
evaluations. This combination, employing worth-category Plan 1 for parts, provides the fol- 
lowing categories. 

1A, 1B, 1C, 1D* 

2A, 2B, 2C, 2D* 

3A, 3B, 3C, 3D* 

4A, 4B, 4C, 4D* 
The numbers refer to the worth of the components; the letters reflect the worth of the part. 
Thus, for example, items falling into Category 1A have the highest valuation. They may be 
described as items where a part shortage cannot be compensated for by on-board manufacture 
or jury-rigging and where failure to maintain the parent component in an operable condition 
would result in the termination of the ship's mission. Items falling into Category 4C have the 
lowest valuation. These may be described as items where a part failure would have no effect 
on the operation of the parent component and where the failure of the component itself would 
have a negligible effect on the accomplishment of the ship's mission. 

For Table 6, the part (based on Plan 1) and component valuations have been combined to 
provide for all categories of material studied, excepting equipage and those items which can- 
not be installed by the ship's force. 

Table 7 was constructed on the assumption that Categories 1A, 1B, 2A, and 2B include 
the high-value items. Items falling into the remaining categories, under this assumption, are 
considered to be of low military value. 





*Military-worth Category D refers to items which cannot be installed by the ship's force. 
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TABLE 6* 
Distribution of All Parts Pertaining to Each Type of 
Component Assigned to the Combined Worth Cate- 
gories Employing Plan 1 














Worth 
Category Mechanical and : 
Electrical Electronics 

1A 1857 0 
1B 434 0 48 
1C 229 0 17 

2A 1498 1123 
2B 227 1666 61 
2C 253 80 49 
3A 3292 2829 89 
3B 415 5921 14 
ures of 3C 523 405 71 
comes 4A 767 1107 293 
vorth 4B 104 1714 176 
2 fol- 4C 253 230 101 




















= 9852 15,075 1410 





*Equipage and items which cannot be installed by the ship's 
force have been excluded. 








TABLE 7 

wa Percentage of Items for Each Type of Component Assigned to 
has High- and to Low-Worth Combined Categories 
facture 
wae Part (percent) 
ive the Worth 
effect P 

Category Mechanical and 
vould Electronics Ordnance 


Electrical 








ined to 


High 40.8 18.5 42.5 
h can- 


Low 59.2 81.5 57.5 




















include 
on, are MB SIMMARY AND CONCLUSIONS 
The purpose of this study was to determine the feasibility of obtaining a relative meas- 
teof the seriousness of spare-part shortages. The study was primarily directed at develop- 
a measure of military worth (essentiality) which can be used as one of the several param- 
“ers in the determination of shipboard stock levels. 


The approach was based upon a consideration of two interdependent factors as influencing 
‘seriousness of part shortages. One of these is mission effect, which measures the effect of 
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component failures on the ship's ability to execute its assigned mission. The other, mainte. 
nance potential, has to do with the effect of end item or part failures on the operability of the 
parent component. Where such failures would render the parent component inoperable, maip. 
tenance potential considers the capability of the ship's force, in the event of a part shortage, 
to maintain the component in a satisfactory operable condition through on-board manufacture 
of the required part, and/or cannibalization, and/or the employment of jury-rigging proce- 
dures. Part shortages are relatively more serious when, in the event of such a shortage, on. 
board manufacture, cannibalization, and jury-rigging are not feasible. Additionally, part 
shortages are serious to the extent that parts are related to components whose failures would 
impair the ship's capacity for fulfilling its mission requirements. 

The development of classifications or categories of spare parts essentiality results 
from the combination of these two elements—mission effect and maintenance potential. Each 
spare part is ultimately assigned to one of these military-worth categories. 

There are two important results of this study. The first is the large amount of agree. 
ment among the independent answers for the component and part questionnaires. In the case 
of the components, substantial agreement among the three participants was obtained for 1132 
(or 92.4 percent) of the total number of components evaluated. 

The second result of interest is the count of the numbers of components and parts 
falling into each worth category. In the case of components, four categories of worth were 
used. Components of the highest worth were assigned to Category 1, those of the lowest worth 
were assigned to Category 4, etc. Approximately 3 percent of the total number of components 
evaluated were assigned to Category 1; 12 percent were assigned to Category 2; 43 percent to 
Category 3; and 42 percent to Category 4. These findings are extremely interesting in their 
marked contrast to the widely held assumption of equal worth for all components at all times 
insofar as allowance-list preparation decisions are concerned. 

Insofar as the military worth of spare parts is concerned, an over-all average of 66 
percent of the parts fall into the lower worth categories. This indicates that the range of tech- 
nical items vitally necessary for a vessel to fulfill its mission requirements is quite limited. 

There are some particular important advantages to the type of approach employed. (ne 
is that it permits an easy handling of worth information from the standpoint of data-processing 
requirements. A second is that the length of time necessary for the participant to answer the 
questionnaires is relatively short. The approach lends itself to implementation on a routine 
basis. This military-worth technique, with slight modifications in the questionnaires, is 
equally applicable to ships other than submarines. 
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SOME PROBLEMS IN QUEUEING WITH DYNAMIC PRIORITIES*t 
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This paper presents some mathematical and numerical results con- 
cerning certain one-machine discrete-time queueing models, subject to 
various queue disciplines. Interest is focused upon systems in which 
service order is governed by ''dynamic priorities,'' under which a cus- 
tomer's relative priority is in effect upgraded withthe passage of wait- 
ing time, so that the proportion of newcomers who can "'step in front of 
him'' becomes smaller and smaller as he waits. 











NTRODUC TION 
Queueing systems in which the place in line taken by an entering customer depends upon 
his "urgency class" and those of other customers already present, and also upon how long the 
latter have been waiting, are of great practical importance but have received little attention 
from theorists. This paper presents some limited results concerning certain single-server, 
discrete-time systems of this kind. Analogous continuous-time models are familiar in two 
extreme cases: the "first-come, first-serve" case, in which customers are served in order of 
urival, and the "'static-priority" case, in which urgency classification alone is the basis for 
sequencing customers (except when ties occur). These cases are considered below, but the 
main purpose of this paper is its treatment of the "dynamic-priority" case proper. In this 
(ase, entering customers are assigned to numerically labeled urgency classes, according to a 
lixed probability distribution; and a newcomer takes precedence over a customer in the queue 
and only if, the former's class is numerically superior to the latter's by an amount at least 
wal to the length of time the latter has already spent waiting. This sort of queue discipline 
tilects decision rules used in a variety of priority-governed systems where management's 
meer with a given customer increases as he waits. Examples occur when production jobs 
‘processed in order of promised delivery dates, in a fancy restaurant's treatment of regular 
Moccasional customers, sometimes in the sequencing of aircraft takeoffs and landings, etc. 
The techniques of this paper can be adapted to a variety of problems, some of which are 
cerned with continuous-time systems. However, the simplest applications and those which 
Mi themselves to the most extensive sampling of numerical results seem to me of widest 
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interest, and the ensuing subject matter is limited accordingly. Also, for the sake of brevity 
and simplicity, tiresome (but essentially elementary) mathematical details are often sup- 
pressed, and appeals are made to the reader's "insight" as well as to his purely mathematica 


judgment. 


MODELS 
The three cases of this paper relate to models which are identical, except for the rules 
used to make certain choices among waiting customers. They are all described within the {o}- 


lowing framework of chance events, which occur in the order listed at each integer "clock time’. 


1. Unless the system is empty, with probability q a service is completed, and the 
customer then in service is discharged. 
. With probability p, p< q, a customer enter the system and is randomly assigned 
an integer "urgency number," j, according to the distribution, Pr(j < n) = TH 
. If the system was previously empty and a customer has just now entered, or if a 
a customer has just now been discharged and the system is not empty, then some 
customer presently in the system goes into service. Three cases are distinguished 
according to which of the following rules is used whenever a choice among cus- 
tomers is possible: . 
First-come, first-serve. Choose the customer with minimum arrival time. 
Dynamic priorities. Choose a customer for whom the difference between urgency 
number and time interval since arrival is a minimum, breaking ties in favor of the 
customer with minimum urgency number among those tied. 
Static priorities. Choose a customer with minimum urgency number, breaking ties 
in favor of the customer with minimum arrival time among those tied. 
(Since two customers cannot arrive simultaneously, each of these rules invariably 
determines a unique choice.) 
The dynamic-priority case "approaches" the first-come, first-serve case as the urgency- 
number distribution shrinks to a single value and "approaches" the static priority case as this 
distribution is spread out further and further. 

The wording above is intended to imply that the server works on at most one customer 
at a time; always works if any customer is in the system; and, having started on one customer, 
finishes before turning to another. Also, the three types of chance events are assumed to be 
mutually independent and independent of past history. 

If p, q, and the time unit are small, then the systems specified above approximate 
continuous-time systems with Poisson input processes and exponentially distributed holding 
times. Several unpublished results for such systems agree in "reasonable" limiting senses 
with the results of this paper. 











NOTATION 
The following notations, including those defined previously, will be used throughout: 
p probability of an arrival, in any given cycle. 
Th probability that a given entering customer will have an urgency number not 
exceeding n. 
q probability of a completion, in any given cycle (if the system is not already 
empty). 
p/q = proportionate "utilization" of the server. 
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= pt, = proportionate "utilization" of the server by customers with urgency num- 





" bers not exceeding n. 
ical 1-4 
“T-. 
pZ 
p/4 
rules (1 - a) 1-p 
fol It follows from earlier assumptions that p, Pa» Z, and A are less than unity. The last three 
t time" notations listed will be given meaning later in the paper. 
In some cases a mnemonic consistency will be maintained by the use of one letter, 
variously "decorated," to refer to related quantities in the three cases; e.g., W is the equi- 
ned librium mean waiting time in the first-come, first-serve case, W: is the equilibrium mean 
waiting time for customers with urgency number j in the dynamic-priority case, and W. is 
- the equilibrium mean waiting time for customers with urgency number j in the static-priority 
” | case. 
lished 
< RESULTS VALID FOR ALL THREE CASES 
The following assertions hold for all of the systems considered in this paper. 
ASSERTION (1): Let Ss). be the equilibrium probability that there are just k customers 
ped inthe system at a "low point''—just after a customer has been discharged, if this occurs in a 
of the given cycle, but just before a new customer enters, if this occurs. Then 
ig ties S. i (1 - A) ak. 
iably ASSERTION (2) Let X(t) be the equilibrium probability that the total service time to 
te required in the future, by all customers present in the system at a given "low point," will 
‘y- teless than or equal to nonnegative integer t. Then 
as this 
Xo(t) = 1 - pzt 1 
stomer 
stomer, ASSERTION (3): The over-all equilibrium mean waiting time, for customers of all 
to be wgency classes taken together, is W. 
Assertion (1) is established by solving the following "equations of detailed balance," 
ate long with the equation us, = 1: 
Iding 
nses 
p (1 - q) Sg = q(1-p) sy 
[p (1 - q) + q (1 - p)] 5, = p(1-q)s,_,+@(1-p)s,\4 
= k= 1,8 3,... 
ot 


or examples involving essentially the same technique, see W. Feller, Probability Theory and 


ty = (Wiley, 1950), or P. M. Morse, Queues, Inventories and Maintenance (Wiley, 
8), 








| Assertion (2) can be thrashed out of the observation that the random variable involved 
Sisum of mutually independent, geometrically distributed variables, each with mean 1/q, 
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and the number summed being k with probability S,- However, it is algebraically easier to 
analyze a pertinent stochastic process. Let x, (t) be the equilibrium probability that, after the 
customers who are in the system at a given "low point" have received just t units of service 
time in all, i of them will still be in the system. This is consistent with the definition of Xo(t) 
in (2); and (1) asserts that x, (0) = (1 - A) A’. Elementary probability considerations yield the 
following relations, for t= 1, 2,3,...: 


X(t) = X(t - 1) + q x, (t - 1) 


x, (t) = (1- aq) x,(t- 1) + ax,,, (t- 1), toi 2a c.. 


(See Feller, op. cit., and Morse, op. cit. , for similar examples.) It is a straightforward matter 
of induction to show that x,(t) = (1-A) A’ z', for t= 0,1,2,..., and i= 1,2,3,.... Finally, 


co 


xp(t) = 1- ». x,(t) = 1- pzt*t, 
i=] 


as required. 

As a step toward establishing (3), it is easy to show directly that W is the mean of the 
distribution of (2). The reader should be able to convince himself that the over-all mean wait- 
ing time must be equal to the mean total service time for customers in the system at the point 
in the cycle when arrivals occur, regardless of the way in which customers are selected from 
the queue, provided that the manner of choice is independent of the amount of time to be 
required for service (under the conditions imposed on all of the systems treated in this paper). 


FIRST-COME, FIRST-SERVE AND STATIC-PRIORITY "EXTREMES" 

The following assertion, giving the equilibrium waiting time distribution in the first- 
come, first-serve case, is an immediate consequence of (2). 

ASSERTION (4): In the first-come, first-serve case, the equilibrium waiting-time 
distribution, over nonnegative integers t, is given by 


Pr (wait <t) = t-sf 1. 


The static-priority case lends itself to a treatment precisely parallel to that used by 
A. Cobham (Operations Research 2, 1, 1954) to solve a variety of continuous-time static- 
priority problems. The proof of the following assertion will not be reproduced here, since it 
can so easily be reconstructed by reference to Cobham's article (and several other discrete- 
time results can be established in the same way). 

ASSERTION (5): In the static-priority case, the equilibrium mean waiting time, Wj 
for customers with urgency number j, is given by 





» 7 1-p 
Wo (1 - By) A - A)” 
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DYNAMIC PRIORITIES: BOUNDS ON EQUILIBRIUM MEAN WAITING TIME 

The bounds given in (6) and (8), below, indicate how the dynamic-priority case relates 
to the "extreme" cases above, and (7) will be computationally useful later. 

ASSERTION (6): In the dynamic-priority case, the equilibrium mean waiting time, W; a 
for customers with urgency number j, is bounded as follows: 


Wj 


(6a) + 


= 


(6B) 


¥; 
(6c) —<1+(1-Z) = p, 
Ww icj 3 


ASSERTION (7): In the dynamic-priority case, suppose there is a maximum urgency 
number, J. Let xo(t) be the equilibrium probability that a customer with urgency number J 
will wait t units of time or less for service, t a nonnegative integer. Then 


' 


W,=W+ = p,[i-x,(J -1- i)]. 
J i<J j [Xp 


ASSERTION (8): In the dynamic-priority case, suppose there is a maximum urgency 
number, J. Then Wy is bounded as follows: 


=1+(1-2) 2 p20 
i<J 


Assertion (6A) is established by noting that if the system under consideration is 
replaced by one which is identical, except that all urgency numbers less than j are increased 
to j, and if this system is replaced in turn by a static-priority system with the same p,q, and 
urgency-number distribution, then neither change can increase the mean waiting time for cus- 
tomers who "originally" had urgency number j. The right-hand side of (6A) is the mean wait- 
ing time, obtained from (5), for customers with urgency number j ‘in the twice-modified system. 


The proofs of the other assertions above are closely related. It is useful first to observe 
that 


co in’ 
w-u+) am)” p., 


t=1 =j-t 


Where U is the equilibrium expected delay to customers with urgency number j, due only to 
customers already in the system when they enter, and dw, is the equilibrium probability that 
such a customer will wait exactly t units of time before being served. Now U = W, regardless 
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of j; and since no "inside" summation can exceed tpiiy , the double sum cannot exceed 
Pi-4 W; . Substitution and rearrangement lead to (6B). 

Reversing the order of summation in the expression for W; displayed above results in 
the following: 


Relation (6C) follows from simple manipulations after substitutions based upon the observations 
that U =< W and that no "inside" summation can exceed the probability of a nonzero wait—which 
in turn cannot exceed the probability of a nonempty system, which is A, according to (1). 
Assertion (7) is no more than the specialization of the last-displayed expression for j = J, 
taking into account that in this instance U= W. To obtain (8), note that since a customer with 
urgency number J must wait at least for all customers in the system when he enters, his 
probability of a wait not exceeding t units of time cannot exceed the probability that all the cus. 
tomers in the system at a "low point" will require a total of no more than t units of time for 
service. Then, from (2), Xo(J -1-i)=1- pz} , and (8) follows, after simple manipulations, 
from the obvious substitution in (7). 

Remarks: The following are reformulations of (6A) and (6B), by use of the notation 
defined in (5): 


1-p 


= 


a: 


These expressions, compared with the original versions of (6A) and (6B), indicate how, ina 
very rough sense, the dynamic priority case falls "between" the two "extreme" cases. 

The detailed proofs of (6) and (8) actually make it clear that, if certain degenerate 
cases are excluded, then the inequalities given can all be rewritten as strict inequalities. 

A lower bound for nonmaximum urgency numbers can be obtained from (8) in a manner 
related to the proof of (6A). The general expression will not be given. The technique required, 
which can easily be applied directly to individual problems, is to compare the system at hand 
with one which is identical, except that customers with urgency numbers greater than j are 
not served. In this new system, j is the maximum urgency number, so (8) applies. Also, cus- 
tomers with urgency number j will spend no more time waiting in the modified system than in 
the original one, so the lower bound obtained also holds for the original system. The bounds 
so obtained, incidentally, always improve upon (6A) for some nonempty range of urgency num- 
bers at the upper end of the scale. 


Figure 1 illustrates the strengths, such as they are, and the more apparent weakness¢ 
of (6) and (8). 


WAITING-TIME DISTRIBUTIONS 

We now present a procedure for calculating successive waiting-time probabilities in 
the dynamic-priority case, conditional upon the state of affairs at the time of a customer's 
entry into the system; point out how this procedure can be used to determine the equilibrium 
waiting-time distribution and its mean for customers in a least-urgent class (maximum 
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possible urgency number); suggest how it can be used MEAN WAIT 
to obtain improved bounds on equilibrium mean wait- 
ing time for nonmaximal urgency numbers; and indi- 
cate how the procedure can be modified slightly for 
application to the static-priority case. 

Consider a customer with urgency number j. 
Let Vy(t) be the probability that there will be exactly 
k customers present inthe system whose service will 
precede service to the customer under consideration, 
when t units of time have passed since his entry, k, 
t=0,1,2,... . For each t, let y(t) be the infinite 
column-vector whose successive components are 
pt), yy), yo(t),... . Define the infinite matrix, 





























Figure 1 - Bounds on equilibrium 
R(t), by mean waiting time, W; , by urgency 
class, for the dynamic priority 
example with p = 0.4, q=0.5, and 
the customers equally divided among 
urgency numbers 0, 2,4, and 6, The 
circled dot on j = 6 gives the actual 
mean waiting time for customers 
with this (maximum) urgency num- 
ber, manually calculated by use of 
(7) and the procedure in the section 
on Waiting-Time Distributions. The 
broken line shows W, for the static 
priority example with the same p, 
q, and urgency-number distribution. 
It should be remarked that (6B) is 
not always weaker than (6C) and 
tends to be stronger for (relatively) 
small urgency numbers when the 
= pz, t (1-q), entire range of urgency numbers is 
J- large. A simulation has shown that, 
with errors probably not exceeding 
Q=4 (1 - PT; 4) . 0.2, Wo = 2.6, W2 = 3.5, W4 = 4.4, 
We = 5.4. 











and blank spaces represent zeros. The following relationships can be obtained through 
elementary probability considerations: 


y(t) = R(t) y(t-1), t=1,2,3,... 


knesses These relationships permit the calculation of arbitrary y(t), conditional upon y(0), the vector 
whose components are the probabilities that the customer under consideration will enter the 
system to find present in it 0, 1, 2,... customers whose service will precede his. The prob- 
ibility that the designated customer will wait t units of time or less for service is precisely 

ies in Hi), 80 the recursion permits the determination of as many elements of the waiting time dis- 

er's tribution as may be desired—again conditional upon y(0). 

ibrium 


m 
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Least-Urgency Class Under Equilibrium 

If there is a maximum possible urgency number, J, then the y(0) needed to Carry 
through the calculations to determine equilibrium waiting-time probabilities, for customers 
with urgency number J, is given by (1): 


y,(0) = (1- a) AX, k=0,1,2,.., 


It is not difficult to see that if the cumulative waiting-time probabilities, Yo(t) » are needed for 
t=0,1,..., T-1, they can be obtained by a calculation involving about 3T” arithmetic oper- 
ations, and that obvious efficiencies are possible in computations upon a group of problems 
which have the same R(t) in common for the first "few" values of t. 

The Yo(t) calculated for the maximum-urgency-number class under equilibrium are the 
X(t) called for by (7), as bases for determining the equilibrium mean waiting time, W:. If 
there is a minimum urgency number as well as a maximum one, then only finitely many of the 
Yo(t) are needed to get a precise answer from (7). Specifically, if the range of urgency num- 
bers is T, then the Yo(t) = x(t) are needed for t=0,1,..., T-1. 


Bounds 

The opportunity to calculate equilibrium mean waiting times for a least-urgent class 
can be used to obtain bounds applying to other classes, which often will improve upon the 
results given in the section on Dynamic Priorities, especially for urgency numbers j with p, 
close to p. A lower bound on W. is the equilibrium mean waiting time for customers with 
urgency number j in a system exactly like the original one, except that customers with urgency 
numbers greater than j are not served. An upper bound on W, is the equilibrium mean waiting 
time for customers with urgency number j in a system exactly like the original one, except 
that all urgency numbers greater than j are reduced to j. 





Static-Priority Modification 

The procedure of this section is modified to apply to the static-priority case by replac- 
ing R(t) in the recursion relationships by R(1). The equilibrium waiting-time distribution for 
a least-urgent class can again be calculated, step by step; but (5), of course, provides a much 
more economical way to get the mean value. 


TWO URGENCY CLASSES 

The foregoing results are more complete in their implications concerning the special 
instance where there are just two urgency classes than for general urgency-number distribu- 
tions, and the greater simplicity of the situation makes it possible to obtain some additional 
results. 

No generality is lost by assuming that the two urgency numbers are 0 and J, where J 
is a positive integer in the dynamic priority case proper—but J =0 and J= o will represent 
the first-come, first-serve and the dynamic-priority cases, respectively. The following nota 
tions will supplement those defined earlier: 

Dp 9 * P- Po = proportionate "utilization" of the server by customers with urgency 

number J. 
w(J) = equilibrium mean waiting time for "preferred" customers (with urgency number 
0) in a system where nonpreferred customers have urgency number J. 
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Q(J) = equilibrium mean waiting time for nonpreferred customers in a system where 
their urgency number is J. 

The redundancy of this notation for J = 0 will cause no confusion. The quantities w(J) and 
QJ), of course, depend upon a set of parameters (e.g., q, p, and Po): as well as on J. 

Assertion (9), below, is a straightforward specialization of material from the section 
on Waiting-Time Distributions. The portion of (10) concerned with Q is the specialization of 
(1), and the portion concerning w follows in the light of (3). 

ASSERTION (9): Let 
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x'(0) Column-vector transpose of (1- A) (1, A, A“,...), 


x'(t) = R(<J) x'(t-1), 
x'(t) = R(>J) x'(t-1), t=J+1, J+2, J+3,.... 
Then X(t) , the first component of x'(t), is the equilibrium probability that a nonpreferred 
‘ustomer will wait t units of time or less for service, ina system where such customers have 
gency number J. 
ASSERTION (10): Let x(t) be as in (9); and let X = Zxp(t), summed over 
201,..., J-1. Then 


w(J) = W-fp (J-X) and QU) = W+ py (J -X). 


The appendix gives numerical values of the functions w and 2. We continue here with 
‘ue general information concerning their character. 
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ASSERTION (11): The function w is nonincreasing and convex downward; w(0) = W; 
and as J approaches infinity, w(J)/W approaches 


Po 
—— « 1 - ———-. 
W 1- Pp 


Further, for all J, w is bounded as follows: 


Po [1 - (Z - 11) Z+ 1) ] 





1 - 
1- 1% 


If Po > 0, then w is strictly decreasing and convex, and for 1 < J < the bounds above can 
be written as strict inequalities. 

ASSERTION (12): The function 2 is nondecreasing and concave downward; 2(0) = W; 
and as J approaches infinity, Q(J)/W approaches 


Further, for all J, 2 is bounded as follows: 


Po [1 - (Z - T)Z + 1)?] 2(3) 


= 21+ pj(l-Z 


1 + 
1- 1% WwW 





J). 


If Po > 0, then Q is strictly increasing and concave, and for 1< J < , the bounds above can 
be written as strict inequalities. 

Assertion (11) follows from (12) in the light of (3), so the ensuing (sketched) proof will 
treat only (12). The functional values for J = 0 and J = o are from (4) and (5). The lower 
bound is the direct specialization of (8). The monotonicity statements follow from (10), in view 
of the fact that 0 <xp(t) = 1, with no possibility for equality unless p= 0. The concavity 
statements are based upon successive refinements of the specialization of (6C), each of which 
is an upper bound on Q(J) for J=J 0» expressed in terms AT 0); the details will be omitted. 
The statement concerning the behavior of 2 as J approaches infinity can be proved by a 
moderate extension of the proof of the appropriately specialized version of (6B); again, details 
will be omitted. 

It remains to establish the upper bound given in (12). Let e(n,J) be the expected waiting 
time for a nonpreferred customer who enters a system in which such customers have urgency 
number J and finds n customers already present. An analysis based upon the procedure given 
in the previous section (or a direct approach of a type used by Feller, op. cit., in dealing with 


certain gambler's ruin problems) leads to the following relationships, in which P and Q are 
as in (9): 


e(1,J) = 1 + (1-P-Q) e(1,J-1) + Pe(2,J-1), 
e(n,J) = 1 + Qe(n-1,J-1) + (1-P-Q) e(n,J-1) + P e(n+1,J-1) 


n= 2, 3, 4,-+° 
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ithe nth of these relationships is multiplied by (1 - A) A", they are summed, and the sub- 
stitution 2(J) = Z (1 - A) A™ e(n,J) is made throughout; then the result (after elementary 
algebraic manipulations) is as follows: 


Q(J) = A + (Z - t9Z + m9) Q(J-1) - Po(1- 4) (1- Z) e(1,J - 1) 


S03, 2, 5.4. 


Astraightforward induction, based upon the fact that 2(0) = W and the observation that 
e(1,J-1) =1/q, leads to the required inequality. If Po > 0 and J = 2, this inequality is 
strictly observed, since then E(1,J-1) >1/q. (The general technique used above can be 
applied to systems with arbitrary urgency-number distributions, to obtain a computing tech- 
nique roughly equivalent in cost and output to the procedure in the section on Waiting-Time 











can 

Distributions. However, I have been unable to find any interesting general conclusions by use 
_ W: of this approach which are not more easily reached otherwise, except those presented above.) 

Remarks: Note that for J=0 and J =1, the upper and lower bounds of (11) and (12) 
are equal, whence «(0) and 2(0), already known to equal W, are precisely specified, and so 
are w(1) = W - PA and 2(1) = W+ PoA. 
The following bounds, which are weaker consequences of the "new" bounds of (11) and 
(12) may be of interest: 
(1-7 0) | 
oJ) ie Po 1-Z 
WwW 1 = 19 
ve can Py E . a "o7] 
= 
f will 1-1 
er 
— For Z close to unity (which must be the case if a reasonably close approximation to a continuous- 
ty time model is concerned), the weakening is slight. 
‘hich The bounds obtainable by specializing (6A) and (6B), along with (3), are implicit in the 
tied lirst sentences of (11) and (12). The following bounds derive from (6C): 
;, w(J) ~ 
etails iw 1 - Po (1 - Z)J 
waiting 
QUI 

e given WwW 
with 


These are invariably weaker than the corresponding "new" bounds of (11) and (12) but may be 


{some interest because of their relative simplicity. In particular, they make it immediately 
Clear that 


are 


Q(3) - wJ) < pA. 


This has a significant interpretation when urgency numbers are considered to be "scheduled 
ns waiting times," so that the dynamic-priority discipline calls for choosing a customer to mini- 
lize the following quantity: 


562850 O - 60-5 
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"scheduled starting time" = arrival time 
ASYMPTOTE 


+ "scheduled waiting time." 
QJ) 





THTTT TT eee TTT 


If such a situation fits the model of the present 
section, then the above inequality implies that 
inevitably the difference between actual mean 
waiting times for two classes of customers vill 
be less than the difference between "'scheduled 
waiting times." 
Figure 2 illustrates (11) and (12) fora 
ASYMPTOTE numerical example, and includes the actual values 
‘ of w and 2, calculated by use of (9) and (10). 
For convenience, these functions are graphed as 
Figure 2 - Two-urgency class case: w curves, although they are really defined only for 
and 2 as functions of J, for p = 0.4, integer values of J. (It should be clear that the 
Reh gn tig ee ae a admission of noninteger values would not add to 
actual values of w(J) and Q(J) are typi- the range of queue disciplines represented, 
wet Bh ce: amefoneernptcn Mite ee ee although it would greatly complicate the notation, 
prepared. expression of assertions, etc.) 











APPENDIX 


This appendix summarizes computations based upon the final section above. The main 
results are in Table A3, whose arrangement was suggested by the fact that the number 


_W- 03) 2Q)-W 
~ W-wle) 26) -W 





is a natural specifier of the relationship of a given dynamic-priority system to the most closely 
related first-come, first-serve and static-priority systems; plus the (computationally empirical) 
observation that if 6 is considered as a function of p, 9, J log (1/Z), and q, then variations 
in q through the range of 0.5 downward seem to have little effect on 0. 
Tables Al and A2 facilitate specific numerical calculations in which Table A3 is used. 
It was found computationally that the quantity, 


1- (1- p) 2) - p(Z-19Z+ 1), 


provides close approximations to @ for 7 0 < 0.8. This quantity is a "natural" compromise 
between the upper and lower bounds on @ implied by (11) and (12). 





SOME PROBLEMS IN QUEUEING WITH DYNAMIC PRIORITIES 


TABLE Al 
Values of W for various combinations of p and q 





q 





0.02 0.01 








12.25 24.75 
32.67 66.00 
73.50 148.50 
196.00 396.00 
441.00 891.00 
931.00 | 1881.00 
2401.00 | 4851.00 





























TABLE A2 
Values of 1/[log (1/Z)] for various combinations of p and q 





q 





0.02 0.01 





61.75 124.25 
82.17 165.50 
123.00 248.00 
245.50 495.49 
490.50 990.49 
980.49 | 1980.40 
2450.32 | 4949.84 
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CONFIDENCE INTERVALS OF PREASSIGNED LENGTH FOR 
QUANTILES OF UNIMODAL POPULATIONS** 


Lionel Weiss 


Cornell University 





For the case where the only information we have about a population is 
that it has a unimodal density function, a two-sample method is given 
for constructing a confidence interval of preassigned length for any 
quantile of the population. 











INTRODUCTION 

A common problem is that of estimating the qth quantile of a distribution about which 
nothing is known except that it is continuous. The usual procedure for constructing a confi- 
dence interval for the quantile is based on a sample of predetermined size n and uses as the 
confidence interval the interval between the jth smallest value in the sample and the je 
smallest value in the sample, where i and j are chosen so that 


1 





q., i i-i-1 : 
| r*(s-r))""*(1-s)"! ar ds=B, 
0 


n! 

(i- 1)! G-i-1)! Mm-j)! J, 
where 8 is the desired confidence coefficient. The length of the confidence interval given by 
this usual procedure is a random variable, and when the length is large the confidence interval 
isof limited use. In this paper, a two-sample procedure for finding a confidence interval of 
preassigned length A and confidence coefficient of at least 8 for the q~ quantile of a unimodal 
population is described. 

A probability density function f(x) is called "unimodal" if there is a value Xo such that 
ix) is a nondecreasing function for x < Xo and a nonincreasing function for x= Xo- Xq is 


then called a "'mode."" There may be more than one mode, in which case the set of modes forms 
a interval. 


y 
| f(x) dx = q has a unique solution in y, this value of y being called the ngth quantile" of 
“0 
the distribution. 


It is clear that if f(x) is unimodal, then, for any q in the open interval (0,1), the equation 


Suppose xX, Xo » ++» are independent, identically distributed, random variables and all 
that is known about their common distribution is that it has a unimodal probability density 
tere 


‘Manuscript received July 8, 1959. 
Research supported by the Office of Naval Research. 
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function. If we are given positive numbers A, 8, where 8 <1, it is impossible to construct 
a confidence interval of confidence coefficient of at least 8 and of preassigned length A for the 
qth quantile, if the sample size n is also preassigned. (For if it were possible, it would be 
also possible as a special case to estimate the median of a normal distribution with unknown 
mean and standard deviation, and the result of [1] shows this to be impossible.) In this paper 
we describe a two-sample procedure for constructing such a confidence interval, the size of 


the second sample depending on the observations in the first sample, as in Stein's procedure 
for the normal case [2]. 


DESCRIPTION OF THE CONFIDENCE INTERVAL 

B, 4, q are given fixed positive numbers throughout the discussion, 8 and q both 
being below unity. Q denotes the qth quantile of the unknown unimodal distribution. F(x), 
f(x) denote respectively the cumulative distribution function and the probability density function 
of the distribution. If we take a sample of size n, where nq is not an integer, the qth sample 
quantile is defined as the sample value with exactly [nq] observations below it, where [nq] is 
the largest integer not greater than nq. Since the distribution is continuous, we can ignore 
the possibility of ties among the observations. Z will denote the qth sample quantile. 

If a, y are any two values, both in the open interval (0,1), we define’ N(a,y) as the 
smallest positive integer n satisfying the inequality 





min(1,q+7) 
nl i ylna] (1 - yy2-[na]}-1 dy = a. 


(Inq)! (@-[na]-1)! J ax(0,q-7) 


Such a value N(a ,y) can always be found, since the integral approaches unity as n increases. 
We construct our confidence interval as follows. We choose quantities a, w, r in the 

open interval (0,1) such that aw=8 and r> max(q,1-q). We take afirst sample of m 

observations, where m is the smallest positive integer that satisfies the inequality 


1- [mr™-1 - (m-1) r™] =w. We define L, U as the smallest and largest of the m observa- 
tions, respectively. Then we set: 


U-L 2 U-L 2 


Y = [min r-q, r-(1-q), 


r-(i-d A r-a ;| 


Then we take a second sample of size n, where n is the smallest integer which is greater 
than N(a,y) such that nq is not an integer. Denote by Z the qth sample quantile of the sec- 
ond sample. Then our confidence interval for Q is [Z - (4/2), Z + (A/2)]. 


PROOF THAT THE INTERVAL HAS THE DESIRED PROPERTIES 

First we note that if a, y are any two values in the open interval (0,1), and if 
q-F[Q-(A/2)] = y and F[Q+ (4/2)] - q = y, then if we take a sample of size n = N(a,7) 
such that nq is not an integer, P(|Z - Q| = A/2) =a, where Z is the qth sample quantile. 
To see this, we can write P(|Z - Q| = A/2) as 


Q+(4/2) 


n! a 
({nq])! (n- [nq] - 1)! ai [F(@)} Pa) fa - r(ypPo lal? ¢¢z) az. 








); 
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Making the change of variable y = F(z), and using the assumptions about F(y), we see that if 
n2N(a,y) then the integral is at least equal to a. 

Next we note that if we take a sample of m observations from the population, and denote 
the smallest and largest observation by L, U, respectively, then P[F(U) - F(L) =r] = 
mae (mr™~1 - (m-1) rl , for any r in the open interval (0,1). For any given r in the open 
interval (0,1), 1 - [mr™- - (m-1) r™) approaches unity as m increases. We take r as 
some fixed quantity above max(q, 1-q) and below one, and then choose m to make 
1- (mr™~ - (m- 1) r™] =w, where w is some given value in the open interval (0,1) to be 
specified more exactly. If F(U) - F(L) =r, then it follows that L <Q <U; for if Q < L, 
then F(L) 2q, and F(U) - F(L) <1-q <r; while if Q>U, then F(U) =q, and 
FU) - F(L)=q<r. 

Defining y as before, we now show that if F(U) - F(L) =r, then F [Q+ (A/2)]-q2y 
and q- F[(Q- (A/2)] =y. Clearly, if F(U) - F(L) =r, then F(Q) - F(L) = r-(1-q), and 
FU) - FQ) =r-q. First we examine q - F[Q - (A/2) ]. If L =Q- (A/2), then 
q- F[Q - (4/2)] =q- F(L)= r-(1-q). If L <Q- (A/2), there is an ¥ between Q - (A/2) 
and Q such that 


f(x) = 





A 
a-F(a-3] 
=m, 
2 


Then either (a) f(x) = f(x) for all x = Q - (A/2), or (b) f(x) < {(x) for all x=Q. 
In case (a), 


A A 


FQ) - F(L) =q- FQ - )s f(x) (2 oe L)<q-F(@-2)+ 16 (u - L->) ; 


or 





*) stab a 
2/ U-L 2° 


F(U) - F(Q) =f(%) (U - Q) = f&) (U - L), 
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Thus we have proved that if F(U) - F(L) =r, then 


or r-(1-q) A r-qg A 
a-F(@-3)= min|r--a), os 2" z=4 31. 


Next we examine F[Q+(A/2)]-q. If U =Q+(A/2), then F[Q + (A/2)]- q = F(U) - 
q=r-q. If U>Q+(A/2), there is an x' between Q and Q+ (A/2) such that 


F(@+>)-4 


{(x') < A 
2 





Then either (c) f(x) < f(x') for all x = Q, or (d) f(x) = f(x') for all x => Q+ (A/2). 
In case (c), 


q - F(L) = (Q - L) f(x') = (U - L) f(&'), 





In case (d), 


FU) - F(@) =F(@+>)- a+tt) (u-@ ->)s F(@+3)- a+tte) (U-L 2), 


or 





F(q+4 
Q + 2 
Thus we have proved that if F(U) - F(L) =r, then 


*3;° 3* ~~ 8s £* Bem 
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To summarize the preceding paragraph, we have proved that if F(U) - F(L) =r > 
max(q, 1-@), then both F[Q + (A/2)]- q and q - F[Q - (4/2) ] are not less than 





' r-(l-qd A r-qaa 
y=min| rg, 7-(-9), e « %, 2’ F=9 2) > 


Now we can show that the confidence interval [Z - (A/2), Z + (4/2)] described earlier has 
confidence coefficient at least equal to 8. Choosing w,a@ so that wa = 8, we have 

p[z - (4/2) = Q= Z+ (A/2)]= P[Q - (A/2) = Z = Q+ (A/2)] = P[F(U) - F(L) = r] P[Q - 
(4/2) < Z=Q+ (4/2) | F(U) - F(L) =r] =wa = B. 


GENERALIZATIONS 

It is clear that various other combinations of order statistics from the first sample may 
be used instead of the smallest and largest values in order to set the size of the second sample. 
Also, it is easily seen that several quantiles of a unimodal distribution may be estimated 
simultaneously by confidence intervals of preassigned length, the simultaneous confidence 
coefficient being at least equal to 8. All that is necessary is an obvious modification in the 
formula for the size of the second sample. 


APPROXIMATE DISTRIBUTION OF THE SAMPLE SIZE 
For the sake of definiteness, let us assume that q = 1/2, so that q=1-q. Then 


r - (1-q) 4 


y= min|r-(-a, U-L 2 


Define 6(€) by the equation 


ws 1 
[ (1/N27) exp - —t 
6(e) 


Using the fact that 
a [ex an at =) (") o} ep) 
(x-1)! (n-x)! 0 j=0 j , 


4sum of binomial probabilities, and using the normal approximation to the binomial distribu- 
ton, we find that for small values of y , 


1- 1 
Nia ,y) = al = 9” 5 (1 -a)| : 
r 


Tus if U-L > A/2, 





4(v - 1)? q(1-4) a= (1-a)| 
; 


N(a,y) = 6 
a? [r-(1-@) 
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This means that for certain distribution functions F(x), the expected sample size required by 
our procedure is infinite. 

There are strong grounds for suspecting that any two-sample procedure for this 
estimation problem will have infinite expected sample size under some distributions. For it js 
well known that in the case where the form of the distribution is known, but the location and 
scale parameters are unknown, then in order to construct a confidence interval of predeter- 
mined length for the location parameter, the sample size must be proportional to the square of 
an estimate of the scale parameter. Thus, in the general case, the size of the second sample 
must be essentially homogeneous of the second degree in the differences between the order 
statistics of the first sample. But there are distributions for which the square of differences 
between order statistics has infinite expectation. 

There is some freedom in choosing a, w, and r, and we might try to choose them to 
give relatively high probabilities to small values of total sample size. Of course, the best 
choice of these quantities depends on the true but unknown distribution. Perhaps a three-sample 
procedure would be worth studying: the first sample a small one to help in deciding on the 
values of a, w, and r to be used in the second sample, and the size of the third sample being 
set according to the values observed in the second sample. 
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ON ASSURING SAFETY IN DESTRUCTIVE TESTING*t 
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INTRODUCTION 
This report deals with the use of items which are very dangerous if defective; it is thus 





desirable to have no defectives, at least with high probability. The situation is special in two 
other ways: (a) if the items are used at all, it is known that exactly N will be used, and (b) 
testing of an item is destructive. 

The following procedure is considered for deciding whether or not to use the items. A 
surveillance sample of size N, is tested; if none are defective, the items (exactly N) are used; 
otherwise, none are used. 

Let A(p) be the probability that no defective items are used, where p is the probability 
of any particular item being defective. The event "no defective items are used" can occur 
either because the items are not used at all or because there were no defective items among 
the N used. 


Let A = A(p). It is shown that A depends only on the value of N,/N ( = 86); 


more precisely 


B 
- B 1 
(1) ke1-(7*)) vet 


A is a measure of the guarantee that we have of the use of no defectives. A computation by use 
of Eq. (1) shows that in order for A to be 0.99 (fairly good protection), 8 must be chosen to be 
35, which represents a forbidding amount of testing—to use one item, 35 must be tested. Table 1 
gives several computed values. 

The moral of these computations seems to be that the type of perfection desired cannot 
be attained by the kind of dichotomous testing that has been introduced here. Faith in the safety 
of the items must always be generated in some other way, a way which may involve a detailed 
physical model. 





FORMAL DEVELOPMENT 
The following is a list of symbols used along with descriptive or formal definitions. 


*Manuscript received March 18, 1960. 
tResearch supported by the Office of Naval Research under Task Nonr-266(55). 
}The author is now with Syracuse University. 
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TABLE 1 





~ ~ 


B A A 








0.05 0.18 
0.1 0.29 
0.2 0.42 
0.3 0.51 
0.4 0.57 
0.5 0.62 
0.6 0.65 
0.7 0.68 
0.8 0.71 
0.9 0.73 
1.0 0.75 
































probability of a defective, 
number of items which are needed in use, 
size of surveillance sample. 


probability of using the N items when a decision procedure is adopted in which 
the surveillance sample is used, 


the number of items actually used, 
the number of defective items among the n used, 
A(p) = Prob (ng = 0). 


We note that n= N with prob L(p) 
= 0 with prob 1 - L(p). 


A(p) = Prob (ny =0/n=0) Prob (n= 0) 


Prob (ng = 0/n=N) Prob (n=N) 


[1 - L@)]+ (1 - p)% Le) 


1 - L(p) [1 - (1 - p)§). 


If it is assumed that L is continuous and L(1) = 0, then 
A(0) = A(1) = 1, and A(p) has a minimum for some p with 
0<p <1. 

Let A = A(p). The situation is summarized in Figure 1, 


SPECIAL CASE 

The special case mentioned in the Introduction is the 
case when the N items are used only when all N, items of 
the surveillance sample are non-defective. Thus, 
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N 
L(p) = (1- p) § 


N 
Ate) = 1- (1-p) § [1-(1-p)Y}. 


infinding p the ordinary methods of calculus apply, and it turns out that 





‘ 7£)" oe 
( 


For large 8 and A close to 1, the following approximations are good: 


‘la. * 


where e is the base of the natural logarithms. 
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PROGRAMMING THE SUPPLY OF A STRATEGIC MATERIAL *t 


Part I. A Nonstochastic Model 


Herman F, Karreman! 


Princeton University 





This paper deals with the various ways in which the requirements of a 
certain strategic material could be met in a period of limited war. 
Besides a small stock of ores to draw on at the beginning of the period, 
high-quality foreign ores have to be imported or low-quality domestic 
ores made suitable for the production of alloys. In the latter case, 
however, the plants for upgrading and processing the ores first have to 
be constructed. The objective is to meet the requirements, assumed 
here to be given for the years of limited war, at minimum cost. The 
cost function contains, besides linear, also quadratic terms, positive 
as well as negative ones. Two recently developed computational meth- 
ods have been used to solve the problem. 











NTRODUC TION 


Some materials are more important for the proper functioning of a country's economy 
than others. Present-day developed economies can hardly do without minerals like iron, copper, 
tickel, aluminum, or manganese, to mention just afew. They become "strategic" materials 
vien a country relies heavily on importation to meet the requirements of those materials. 


Under those circumstances there is no automatic guarantee that they will always be available 
a reasonable cost. 


This applies in particular to the case of "limited" war, by which is understood a situa- 
tion between "cold" war with no overt hostilities and "total" war. The Korean war, with 
ostilities confined to a local area and lasting for several years is an example of the type of 
‘limited" war that is considered here. In such a situation it will be difficult to obtain sufficient 
quantities of these "strategic" materials at reasonable cost by importation only. In fact, 
trices (including transportation costs) paid in the last "limited" war period for these imported 
strategic materials were in many cases about twice the normal ones. 

To protect themselves against such occurrences, countries have started to buy extra 
qantities of these foreign ores and stockpile them in normal times. However, the amounts of 
uoney needed to provide adequate protection this way are enormous. Even the resources of 


Py 

‘Mauscript received January 20, 1959. 

is paper resulted from a study made on behalf of the Office of Defense Mobilization under 
Contract Nonr-1858(02) with the Office of Naval Research. 

lam indebted to Dr. Philip Wolfe for the first and to Dr. J. B. Rosen for the second, of the 
WO Computational methods that have been used to solve the problem. In its formulation I 
venefitted from the criticisms of Professor Oskar Morgenstern and the other members of 
Princeton University’s Econometric Research Program, 
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the United States were strained, despite its having large deposits of many of these essentia] 
materials in its own territory and large amounts of money available for the purchase of extra 
quantities of these foreign ores. The stockpiles of imported ores, built up in the past by the 
U. S. Government, will in many cases be sufficient to provide the industry with only a fraction 
of the extra quantities of strategic materials that it will need in a limited war period. 

But there are, in some instances, deposits of these strategic materials in the country 
itself. On the average, however, the quality of these ores is inferior to that of foreign ores, 
They have first to be upgraded to meet the standards set by the industry for processing into 
alloys or they have to be treated in a special way. In normal times this makes them more 
costly than foreign ores (otherwise they would have been used in the past). But the technology 
of upgrading the poorer ores is steadily improving, and domestic ores might very well become 
profitable in a future period of limited war when the prices of imported minerals will again 
be high. 

Whether or not such a limited war situation will develop out of the present one of cold 
war is, of course, an open question. But in order to be prepared, it has been assumed that at 
some future time such a situation will have to be faced. To be precise, it has been assumed 
here that we are on the eve of a limited war which will last for several years.@ Starting out 
from this assumption, the problem then becomes: given a small stock of ore at the beginning 
of the period, how much has to be imported and how much has to be produced domestically 
each year to meet the requirements of a certain strategic material in the following n years of 
limited war at minimum cost? 

A model has been constructed to solve this problem. The various ways in which the 
requirements can be met, as well as the restrictions imposed on them, have found their for- 
mulation in this model. The unknowns in these expressions represent the activities that can be 
performed at various stages of production. Their values are the quantities that have to be 
imported, produced domestically, etc.; they have to be determined by the analysis. The prob- 
lem is then to find that combination of quantities that meets the requirements while staying 
within the restrictions and minimizes the over-all cost of the program. 

The prices and transportation costs used in this model for the importation of foreign 
material are derived from past data. The export prices per unit of quantity were found to 
change with the quantities that were exported. It can be expected that this functional relation- 
ship between unit prices and quantities will extend itself into the future. This leads to an 
objective function with a good many quadratic terms besides linear ones. 

As for the upgrading and processing of domestic material, several new processes have 
been developed in the last four to five years. It can be expected that some of them, when oper- 
ated on a large scale, will produce high-quality ores at costs comparable with those of foreign 
ores in times of limited war. These processes have been incorporated in the model and their 
unit costs in large-scale operation assessed. Here, as in the case of foreign ores, the costs 
per unit of quantity also will be functions of the quantities that will have to be produced. 

The latter relationship is, however, quite different from the one that represents the 
cost behavior of foreign material. In the case of domestically produced ores it must be expected 
that the cost per unit will decrease with increases in the scale of production. Purchasing large! 

quantities of foreign ores will, on the other hand, lead to a rise in the cost per unit of imported 





2This is, of course, a rather restrictive assumption. Normally, more than one political devel- 
opment is possible. The reader is referred to [5] for discussion and treatment of more com 
plicated situations. 
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ore. Consequently, the total cost of domestically produced ore is an expression of linear and 
negative quadratic terms, while that of imported ores is an expression of linear and positive 
quadratic terms. The over-all cost of the program, being the sum of these two total cost 
functions, is therefore an expression containing linear and positive as well as negative quad- 
ratic terms. As before, the problem is to find those values of the activities that would make 
the over-all cost of the program a minimum. 

At the time that this problem was formulated there existed no computational technique 
for solving it. A first step which goes a long way toward the ultimate solution was made two 
years ago when a method [2] was developed that could solve the problem if the negative quad- 
ratic terms of the objective function, related to the domestic production of ores, were dropped. 
Hence, the cost expression that was minimized the first time contained linear as well as posi- 
tive quadratic terms for the importation of foreign or2s but only linear terms for the upgrading 
and processing of domestic ores. Decreases in the over-all cost of the program due to 
increases in the scale at which the latter would be upgraded and processed were therefore 
neglected in this first approach. 

This shortcoming was partially overcome by a second computation with a somewhat 
lower linear cost of domestically produced ores while the cost of imported foreign ores was 
kept the same. The difference between the two outcomes showed that a small decrease in the 
wit cost of domestically produced ores would lead to a substantial decrease in the over-all 
cost of the program. It would, moreover, greatly influence the distribution between foreign 
and domestic ores in meeting the requirements. 

This became even more apparent when a third computation was made, based on another 
method [3] just recently developed. As in the first method the problem is first solved without 
taking the negative quadratic terms in the objective function into account. The technique used 
is different, however, from the simplex technique employed in the first method. Once this 
simpler problem is solved, then the negative quadratic terms are brought into play and the 
search for a lower minimum cost continued. The cost function minimized in this third compu- 
tation was exactly the same as the one of the second computation, only it now included the 
negative quadratic terms for the domestic production of ores. The decrease in the over-all 
cost of the program resulting from the latter, as well as the change in the distribution between 
foreign and domestic ores, again turned out to be substantial. 

A detailed analysis of the results, obtained from the first two computations, can be 
found in [1]. A summary of the results of all three computations is given later in this paper, 
where the computational aspects of the problem are dealt with more specifically. 


THE MODEL 


In examining the various forms in which manganese, the material for which this study 
was made, will be required in a period of limited war, one finds that there are essentially 
three rather homogeneous groups of demands for this material. The first group can use only 
product 1 (ferro-manganese), the second group only product 2 (silico-manganese), whereas the 
third group can use both products. Each of the two final products, product 1 as well as product 
2,can be derived in three different ways from intermediate products. Figure 1, which illus- 
trates the technology that underlies the present model, provides a good understanding of what 
follows. As forthe intermediate product, a distinction had to be made between imported mate- 
tial and upgraded domestic ores of the same high quality on the one hand and upgraded domestic 
wes of medium quality on the other hand. The first group of intermediate products can be 





264 H, F. KARREMAN 


processed along conventional lines (processes 3, 4, and 6) into final product, whereas the sec. 
ond group requires special treatment (processes 5, 7, and 8). Another complication is that 
basic material, flowing from source R,,, can be upgraded (process 13) and then further 
processed into product 1 or 2 or, after having been blended with other source material, 
processed into product 2 by one of the two nonconventional processes. There are still further 
complications, as will be seen later, so that it becomes difficult to find that combination of 
quantities to be imported and produced domestically which minimizes the over-all cost of the 
program—a point which acquires even more force when it is recalled that the costs are non- 
linearly related to their quantities. 


REQUIREMENTS 


FINAL PRODUCTS 


THIRD -STAGE OF 
PRODUCTION 














SECOND STAGE OF 
PRODUCTION 


FIRST STAGE OF 
PRODUCTION 
O PRODUCTION 


PRODUCTION AND 
PLANT CONSTRUCTION 
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SOURCES 


A model has been constructed to solve this problem efficiently. The notation used in 
this model for the various parameters, variables, etc., is as follows: 


(a) Parameters in price and cost functions a, Bp 
(b) Predetermined quantities > 
In each individual year 
In a period of years 
(c) Quantities to be determined > 
(d) Technical coefficients 





3In net tons of pure manganese. 
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(e) Indices 
For current years 
For summation of years 
For more than one process 
For individual processes 
The quantities that have to be imported or produced domestically (the x's), the increases in the 
plant capacities (the y's), and the quantities of domestic source-material to be extracted (the 
z's) can be found in Figure 1 along the lines leading away from the circles that represent the 
corresponding activities. 
To start with the demand for the two final products, they are assumed to be known in 
the three quantities ST , 'S) » and &3 for each of the i= 1, 2,....,n years that the program 
is supposed to cover. This is not such a strong assumption as might seem at first sight, since 
only one political situation, the one of limited war, has been assumed to exist during this entire 
period of n years. The quantities that will then be needed every year can be estimated rea- 
sonably well on the basis of similar experiences in the past.* 
However, what is not known at the outset is the extent to which product 1 as well as 
product 2 will contribute to the satisfaction of the bo. ; requirements. Denoting the quantities 
of these products by Xi and Xo i we have: 


) Xi + Xi = Soy 


The quantities Xy i and Xo i have to be determined by the analysis; once determined, they will 
together just satisfy the demands ¢ 2,i° 

The importance of these strategic materials for the rest of the economy made it man- 
datory that sufficient quantities of product 1 as well as of product 2 will be available at all 
times to satisfy the demands for them. In other words, there is no place in this model for 
shortages of final product. On the other hand, technical reasons, e.g., changes in the demand 
specifications, made it less desirable to produce in any one year more final product than will 
strictly be needed that year. No allowance has therefore been made in this model for the pos- 
sibility of building up stockpiles of final products for use in later years.5 Consequently, the 
quantities of final products that have to be produced each year must be equal to the demands 
for it. Stated more precisely: 


= Crit 13 
8, = Sgi t Xo 


vhere the dummy variables f i and Si stand for the quantities of products 1 and 2, respectively, 
that have to be produced in year i. Their values have to be determined by the analysis, just as 
those of X, ; and Xo. 

> ? 


‘Small deviations from these estimated quantities, as well asfrom some of the other constraints 
inthe system, can be taken into account by the technique of chance constrained programming 
described in [4]. 

This may be accomplished when desired (e.g., to extend this analysis to other materials) by 
replacing certain equations in this model with their corresponding inequalities. 
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Looking at Figure 1, one sees that product 1 can be obtained in three different ways 
from intermediate products; so can product 2. Denoting the corresponding quantities of product 
1 to be produced by Xg i> X4j and X5 ,; we have: 

, > > 


f= Xi t Mga t X55 
Substituting this expression for fi into the previous one yields: 
(2) Xgit Xai t X54 = Sait Xi 
Similarly, we find for product 2 (see Figure 1): 
(3) Xgi t+ % i+ %gi = S35 + X25 
The meaning of the expressions (1), (2), and (3) is the following: the quantities Xj and 


Xo i have to be determined so that together they will be equal to Co | i and the quantity of 
penteet 1 that is to be produced must be equal to ¢ 1,i + xy i and that of — 2 equal to 


can +m 

a Vad making use of Equation (1) it is possible to remove one of the two variables Xy i? 
Xo i , for instance x1 ., from the three equations that have been derived so far. However, the 
ain of X1 i have to ” nonnegative to be economically meaningful. Or, what amounts to the 
same thing, the values of Xo i have to be smaller or at most equal to those of ont Hence, 
the first expression becomes: 


(1') Xo = Soi 
After elimination of x, ;, Equation (2) becomes: 
> 
(2") Xo + %gi + Xai + %5 = Sa + So 


and (3) can be written as: 


(3) “Xo t+ Xi t Xi t Xgi = S35 


There are essentially two ways of obtaining the intermediate products used in the third 
production stage, i.e., importation of high-grade foreign material or upgrading low-grade 
domestic source material. In addition, there also will be some high-grade material in the 
stockpile at the beginning of the program. 

To keep the presentation as simple as possible, it has now been assumed that all high- 
grade material, whether imported or produced domestically, will be added to the stockpile S 
and that all high-grade material needed for the manufacture of final product will be taken out 
of S. (This assumption is not too unrealistic and greatly facilitates the presentation.) It will 
then always be possible to meet the requirements if there is at all times sufficient high-grade 
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material in the stockpile to carry-out the (3), (4), and (6) activities at the desired levels.® 
sated more precisely: 


i 


i i 
D1 %3,5%35 +2, %45%4j +2, %65%j = 8, OH 12, 


j=1 j=l j=1 


The technical coefficients Cy j? C4 j? and Cg j? used in this expression, indicate the 
> , > 


net tons (NT) of manganese in high-grade material needed in processes 3, 4, and 6 in year j 

to produce 1 NT of manganese in one of the two forms of final product; they are all larger than 
1. 8; stands here for the quantities of high-grade material present in the stockpile at the 
beginning of the program plus what will be added to it by importation and domestic production: 


i i 


+2. He +2, 


j= 1 j =] 
Combining the last two expressions yields the relation: 
i 


i i i 
a "as “as * pa "e545 * 2 "S34 * 2, *9,j 
]= 


j=1 j=l j=1 


i 


j=1 j=1 j=1 


i i i 
: >. *10,j ~ », 1157 7 X19 5 7 ». X13) = 54 
ft 


The quantities on the left side of the = sign are all to be determined by the analysis, 


8, being the only predetermined quantity in this expression. 

Due consideration also had to be given the fact that the X4, Xg, and X7, as well as the 
%5 Xg, and X15 activities, consume large quantities of electric energy. And the possibility 
could not be excluded that there would not always be sufficient electricity available in a period 
of limited war to carry these activities out at the desired levels. This made it necessary to 
include the following two constraints in the model (see Figure 1): 


6) ge %qy * (egy * Mey * May 


(6 , ‘ < 
“3 *aa * “6 * "955 * 7a * “ee Os * Ms Oe 3% 

ee 

Ot should be noted here that this is a sufficient but not necessary condition for meeting the 
Tequirements since it is also possible (activities (5), (7), and (8), in Figure 1) to obtain final 


products from domestic-source material directly without reducing the quantity of high-grade 
material in the stockpile. 
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The e's in these expressions are technical coefficients like the c's in the previous 
ones; in this case they indicate the quantities of electric energy needed to produce 1 NT of 
manganese in the form of intermediate or final product. "3 and 195 stand for the quantities 
of electric energy that will be available to these processes in year i. The reason why there 
are two relationships here, rather than one, is that the electricity needed for the first and sec. 
ond group of activities has to come from different sources. 

The next feature that had to be taken into account was that the production facilities 
necessary to carry out the X65 and Xg activities and the X9> 13> Xyq> Xy5> Xyg> Xyg, and 
X19 activities did not exist. The construction of these facilities—to the extent that they might 
turn out to be necessary—therefore had to be programmed as well. This led to the inclusion (12) 
of another eight inequalities in the model, one less than the number of activities involved, since 
the X5 and Xg activities will make use of the same facility. It will, however, take more time 
to produce 1 NT of manganese in the form of product 2 in this facility than it will take to pro- 
duce it in the form of product 1. Denoting the ratio between these two time periods in year i (13) 
by ri, we have: 


(10) 


(it) 


X56 + %g i = B53 > 

(14) 
where Ks i stands for the capacity, expressed in net tons of manganese per year, of the facility 

at the beginning of year i. This capacity will be equal to the capacity at the beginning of the 
program, K5 1 plus what will be added to it up to the beginning of year i. Denoting these addi- 

tions by Y5, i? another type of activity, expressed in net tons of manganese per year, we have: made 


i-1 
Ks5i = B51 + 2 V5 j 


J=85 


where 85 denotes the gestation period, the time that elapses between the beginning and the 
completion of a construction, rounded off to the next whole year. 
Combining the last two expressions yields: 


(7) Mga + yy %es - i ¥5,j = * 
J=85 
conta 
As before, the values of x5 i 5 Xg i? and Y5 ,j at the left of the = sign have to be determined, 
whereas the value of Ks 1 is predetermined. 
There are similar expressions for the seven other construction activities (indicated by (15) 
the inner of two concentric circles in Figure 1): 


i-1 
< : (16 
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 addi- To reduce the number of variables of the model the following substitutions have been 
have: made: 


™* 
X15 
1e X16 = 


* * 


X19 = Cy5 (Cg X5 + Cg XQ) = ce Xs + Cg Xe. 


To assure that the program will not call for more domestic-source material than is 
contained in the deposits, the following set of relations had to be included in the model as well: 
ined, 
i 
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i i 
ih “5,5 *5,j * - 35 %*8j = F191 


j=1 j= 
The substitutions that have been made here are (Figure 1): 
"5 “spas * “9 “ae © "a 
“sy * “35 “me 


S17 * “19 *1g * “1g “7 *% = 7%» 


= ' ois ' t as re 
218 = “19 *18 = “18 “7 *7 = “7 *7> 


Zig = yg Xy9 = Cy (C5 x5 + Cg Xg) = C5 x5 + CQ XQ - 


The dummy-variables Z1g° 247> 217 » 249, and Z19 stand for the net tons of manga- 
nese in the form of source material that will have to be derived from the various sources. The 
right-hand side of (16) indicates that this source material will become available in increasing 
quantities as time progresses; such an increase over time could, however, not be expected to 
occur in the quantities of other source materials contained as they are in certain well-explored 
deposits. 

Finally, there are the usual nonnegativity constraints on the x and y variables to 
prevent them from becoming meaningless in an economic sense: 


(19) 
(20) 


Summarizing, we find that the model involves in general 12 variables of production (the 
x's) and 8 of construction (the y's) for each year that is covered by the program. Some of 
these production variables will have zero values in the first year since the capacity to produce 
will not yet be in existence. On the other hand, all construction variables will have zero values 
in the last year since they would only effect the production in the then following year, which 
lies outside the period that the program is supposed to cover. Hence, the number of variables 
involved is slightly less than 20 times n, where n is the number of years covered by the 
program. 

The values of these n variables are bounded three ways, twice from below and once 
from above: from below because (a) all the values have to be nonnegative (relations 19 and 20) 
and (b) the demands have to be satisfied (relations 2 and 3); from above because the values, 
taken separately or jointly, may not surpass certain limits (relations 1 and 4 up to 18). Within 
these bounds those values of the x and y variables have to be found for which the over-all 
cost of the program will be a minimum. 

The latter consist of the cost of importing foreign material and the cost of upgrading 
and processing domestic material. Both types of cost would be linearly related to their quan- 
tities if the corresponding unit prices and unit costs were independent of those quantities. 
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However, there has not been such an independency in the past, especially with respect to the 
prices of foreign materials. Scatter diagrams of export prices against quantities purchased by 
the United States from foreign countries revealed a functional relationship between these prices 
and the corresponding quantities. This is quite natural from an economic point of view since 

it reflects the dominant position of the United States as buyer of large quantities of this mate- 
rial in the foreign markets. A formula for the export price of, for instance, area 9 in year i 
that was found to fit the data fairly well was’: 


Pg i =*%9,i *9i* %9 i-1 *9,i-1 ~ %10,i *10,i 


- %19 i-1 *10,1-1 ~ %11,1%11,1 ~ %11,1-1 *11,1-1 * 991° 


The first two terms of this expression indicate the extent to which this particular 
export price in some year i of the past depended on the quantities purchased from area 9 in 
the same year and in the previous year. The other four x-terms indicate the influence that 
the quantities bought from the other two areas had on the export price of this particular area. 
That influence was a negative one, as the negative signs in front of these terms indicate, and 
rather small compared with the positive influence of the first two terms in this expression. 
(This is an important point, as will be seen later when the computational aspects will be 
discussed.) Still, it was necessary to include the last four x-terms in the formula to obtain a 
reasonably good fit. Similar expressions have been found for the export prices of the other 
two areas. It was assumed that these price-quantity relationships of the past will also hold in 
the future for the years that the program is supposed to cover. 

Besides the export prices, there were also the costs of bringing the foreign materials 
tothe U.S. As distinguished from the export prices, these transportation costs can be 
assumed to be independent of the quantities of foreign material that will have to be purchased. 
They will fluctuate from area to area in the same way as the freight rates of other materials, 
transported in bulk by tramp steamer. The rates that will be charged in the future for the 
transportation of this particular material from each of the exporting areas to the U. S. in the 
eventuality of another limited war have been estimated on the basis of past experience. 

The combined expression for the import price, being the sum of export price and trans- 


portation cost, for 1 NT of manganese in material to be purchased from area 9 in year i, now 
becomes: 


11 


.* 2, Cy Xi t Mind %i-1) *Foi 
k= 


This formula, in which two thirds of the a's have small negative values, is a shorthand expres- 
sion for the previous one, except for the last term. There Boi was a constant of estimation, 
whereas the Bg in this last expression, being the sum of the previous 8 and the estimated 
portation a in year i, will, just as the latter, fluctuate from year to year. These fluc- 
uations will, however, be independent of the quantities of foreign material that will be imported. 


> accra, 


"The @'s and 8's in this expressions are stochastic parameters and therefore of a different 
lature than the technical coefficients of the previous expressions. Both are, however, subject 
o changes over time. 
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The amount of money to be paid by the U. S. for the purchase of X9 j quantities of 
> 
material purchased from area 9 in year i will then be: 


11 
*9,i Poni ~ *9,i 2, yi Mi + ici *k,i-1) + Agi | G=1,2, 


There are similar expressions for the amounts of money that will have to be paid for 
foreign material purchased from the other two areas. Combining all three expressions into 
one, we obtain for the cost of importing foreign material from all three areas in year i: 


11 11 11 
- Xi Py i = Zz Xi * (Oy Mei + Mid %,i-2) * Be 5 
=9 g=9 k=9 

1 il 


4 4 Z, é i i i 0, + Mit *k,i-1 0,0) 


11 
+ pe Bo i %pi 
t=9 


The cost of all foreign material from all three areas in all n years covered by the 
program will then be: 


11 
(21) Xp 


,i Pei 


11 
ye x i Ki 0,1 + Oy i-1 *k,i-1 *0,i) 


n 


> NT 


i= 


This final expression for the cost of importing foreign material contains 3n linear 
and 18n quadratic terms, of which 12n have negative and 6n have positive coefficients. 

The extraction of source material from domestic deposits (processes 16, 18, and 19) 
will take place under decreasing returns to scale, a common feature for this kind of activity. 
This fraction of the total cost function is therefore strictly convex throughout the relevant 
range so that the cost per unit of product rises when the scale at which it is to be produced is 
increased. This leads to the following expression for the unit cost of extracting source mate- 
rial from, for instance, deposit 16 in year i: 
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Pi6,i = %16,1 716,i * 916, » 


where 16 i > 0 and Pigi is the corresponding cost per net ton, which is to be represented 
in the function that is to be minimized. There are similar expressions for the unit costs of 
source material extracted from the other two deposits, Rig and Rig : 

The situation is different, however, with respect to the cost of upgrading these source 
materials. There the unit cost can be expected to decrease when the scale of production is 
increased. For the unit cost of upgrading source material by, for instance, process 12 we 
have therefore: 


Py2.i = 12,4 *12,4 * Pia » 


where 0125 <0. Again, there are similar expressions for the costs of the other upgrading 
processes, "13, 14, and 15. 

In order to keep the number of variables of the model as small as possible, the corre- 
sponding extraction and upgrading activities have been combined (16 with 12, 18 with 14, and 
19 with 15). This led to combined expressions for the corresponding unit costs. For instance, 
let P12 i be the unit cost of domestic material to be extracted by process 16 and upgraded by © 
process 12 in year i. Since it will take C16 NT of manganese in extracted product to produce 
1NT of manganese in upgraded form, we have: 


Pi2.i = P12 + ©16,i Pi6,i 
= 124 *12,5 + 810.4 + ©1631 16,5 216.1 * 216.1) 
= 19 i %12.i1 + °16 4% 16,1 216,1 * F124 + 16,1 916, - 


Since 21617 165 X16 i= °16,i °12,i 12,1 = 12,1 *12,1 95 in (15), this becomes: 


Pyoi = (y0.4 + %16, 4 °16,1 °12,i) *12,1 * F121 * °16,1 916, 


? ' 
= 19 5 %104 + Fiai> 


ware O19 i >0. The negative effect of 19, i, is more than compensated by the Ox. i°16 i 
“M term (2 i? “16 i? and therefore also On being larger than 1), which makes "the 
"i coefficient positive. The same holds for the ay 4,i and O15 i coefficients in the cor- 
responding expression for Py 4, and Pi5 i However, tee the unit cost of material flowing 
{rom source Riz , there is no ‘such cmiaibitin 


Pi3,i = Pi3,i = %13,i *13,i + Fi3i> 


ud here 13 4 remains negative | 
In onlin to arrive at the total cost of upgrading domestic material, these unit costs 
lave to be multiplied by their corresponding quantities and summed up over the years: 
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This is an expression in 4n quadratic terms (n of which are negative) and 4n linear terms, 

However, some of these costs, namely those for which k = 14 and k = 15 have been combined 

with the costs of the third production stage, again with the purpose of reducing the number of 

variables of the model. Consequently, only 2n of the quadratic (n of which are negative) and 

only 2n of the linear terms of this Ey expression are included in the over-all cost function. 
For the cost of the third production stage we obtain in a similar way: 


= “ki Pi i 


1 2 ' 
(Oy i i + Fic i i) 


This expression now includes the cost of extracting and upgrading material that might come 
from the Rig and Rjg deposits. Contrary to expression (22), in which all a's except a3 
had positive values, all the @'s in this last expression have negative values! The reason is 
again that the unit cost of processing material will change inversely as the quantities to be 
processed. This is true for the conventional processes 3, 4, and 6 and in particular for the 
newly developed methods 5 and 8 as well as for 7. Here, considerable economies of scale can 
be expected to appear with increases in the scale of production. This will more than compen- 
sate for the increase in the unit cost of the preceding production stages (18 and 14, 15 and 19) 
as formulated in (22). The result is an expression of 6n quadratic terms, all of which have 
negative coefficients and 6n linear terms. 

Finally, there is the cost of constructing the prodcution facilities: 


n- 
oe be (ey stay Xe i) ° 


i=l k=5 
12, 13, 14, 15, 
16, 18, 19 


As explained before, there is no need for increasing the capacities of the plants in the 
last year (year n). Consequently, the first summation is carried out over n-1 rather than n 
years. Also, since the processes 5 and 8 will make use of the same production facility, only 
one of them had to be included in the second summation. All the a's in this expression, like 
those in (23), will have negative values for the same reason as was given there. Hence, 
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this expression consists of 8(n-1) quadratic terms, all having negative coefficients and 
8(n-1) linear terms. 

The over-all cost of the program covering all n years is merely the sum of the last 
four groups of costs (expressions (21), (22), (23), and (24)). The total number of negative and 
positive quadratic as well as linear terms is summarized in the following table: 


Quadratic Terms 


Negative Positive 
Coeffic. Coeffic. 





(21) E, 12n 6n 3n 
(22) Ey (remaining terms) n n 2n 


(23) E3 6n _ 6n 
(24) E4 8(n - 1) — 8(n - 1) 
27n-8 Tn 19n-8 





In the process of minimizing this cost function, the negative quadratic terms, in partic- 
ular, created some difficulties, as will be explained in the next section where the computational 
aspects are discussed. 

To complete the description of the model, two minor points of a technical nature still 
have to be made. First, it should be noted that the various costs of the program will occur in 
different years and have to be put on a common basis. Therefore, all expenditures have been 
brought forward to the beginning of the first year. 

Second, there is the possibility that the domestic production processes, once operated 
ona large scale, will have become so efficient during the course of the limited-war period 
that they will be able to face the competition of imported materials when times are normal 
again. This has happened many times in the past. In that case, it would not be correct to 
charge the program for the full cost of constructing the domestic production facilities, par- 
ticularly not for those made at the end of the period covered by the program. Therefore, it 
was decided that the program would be charged only for the depreciation of these facilities 
during the years of limited war. To make sure that these charges would not be too low, the 
life span of these production facilities was taken to be fairly short, a common practice in 
cases of newly developed processes. 


COMPUTATIONAL ASPECTS 

The program was supposed to cover a limited war period of 6 years. Consequently, the 
objective function that had to be minimized contained 154 negative and 42 positive quadratic 
terms besides 106 linear terms (see previous listing for n = 6). The computational technique 
that can solve this problem in a straightforward manner still has to be developed. In partic- 
ular, the negative quadratic terms in the objective function make it hard to find the absolute 
minimum of this cost function. 

But even without these negative quadratic terms the problem could not readily be solved 
at the time it was formulated (the beginning of 1957). The then existing computational tech- 


niques for solving quadratic programming problems still had to be further developed to handle 
problems of this type and size. 





H. F. KARREMAN 


A first step going a long way towards the ultimate solution of the problem was made jn 
the spring of 1958, when the simplex technique, which had already proved its usefulness in the 
solution of large linear programming problems, had been further developed [2] so that it 
could also handle quadratic programming problems. However, in order to be applicable, the 
matrix of coefficients of the quadratic terms had to be positive semidefinite. That meant in 
this case that all negative quadratic terms in the expressions (22), (23), and (24) had to be 
dropped to render the problem computable. In fact, all quadratic terms in these three expres- 
sions were eliminated, as well as the n terms with positive coefficients, in order to avoid dis- 
crimination between the various domestic production processes. On the other hand, the matrix 
of coefficients of all quadratic terms in (21) was positive semidefinite, despite the fact that two he 


pro 





computation was made, 





thirds of these coefficients had negative values. The latter, however, were small compared imp 
with the positive coefficients of the other quadratic terms in this expression, as has already neec 
been pointed out earlier. Consequently, the objective function that was minimized in this first units 
computation contained all the quadratic and linear terms of expression (21) but only the linear units 
ones of (22), (23), and (24). It was clear that, by the deletion of the negative quadratic terms 
of these last three expressions, the result would be biased towards importation and away from mill 
domestic production. Nevertheless, it turned out that in the minimum-cost solution more than alloy 
half of the manganese that would be required in excess of the quantity contained in the initial of th 
stockpile should be produced domestically. The actual quantities, expressed in units of 1000 NT lion 
of manganese in the form of alloys, were found in this case to be: shad 
addit 
(a) required in all 6 years 5330 comy 
(b) contained in initial stockpile 1755 rome 
3575 also | 
(c) to be obtained by importation 1557.4 grad 
(d) to be produced domestically 2017.6 
The over-all cost of the program would amount to $1194.6 million, of which $450.1 mil- a 
lion would be for the processing of high- and medium-grade material into alloys. The remain- 07081 
ing $744.5 million consists of $279.0 million worth of high-grade material contained in the initia 
stockpile (evaluated at its shadow price, the meaning of which will be explained later), $222.0 er a 
million for the importation, and $243.5 million for the domestic production of the additional orogr 
amounts of required material [1, p. 38]. less 
A second computation was made to find out, among other things, to what extent the solu- condi 
tion would be affected by lower cost of upgrading domestic material because of improvements large: 
in technology. The method used in this second computation was the same as before.® The tate t 
objective function that was to be minimized contained again all the terms of expression (21) but be wis 
only the linear ones of (22), (23) and (24). Moreover, the unit costs of foreign material, func- 
tionally related to the quantities that would have to be imported, were the same as those used linear 
in the first computation. The unit costs of upgrading domestic material, here again assumed lave » 
to be not functionally related to the upgraded quantities, were set 5 to 10 percent below those 
used in the first computation. The result was not only a sharp decrease in the cost of all high- nestj 
grade material but also a radical shift in the distribution between imported and domestically terms 
Mot be 
8The technique of parametric programming, so useful in solving linear programming problems. Wed j 
had not yet been developed for problems with quadratic objective functions at the time that this tater 
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produced material. The actual quantities, again expressed in units of 1000 NT of manganese in 
the form of alloys, were now found to be: 


(a) required in all 6 years 5330 
(b) contained in initial stockpile 855 
4475 


(c) to be obtained by importation 779.4 
(d) to be produced domestically 3695.6 


in comparing this result with the previous one, we note a marked difference in the shares of 
importation and domestic production in the acquisition of high-grade material that would be 
needed in excess of the initial stockpile. Of the 3575 extra units needed in the first case, 1557.4 
wits = 43.5% would have to be imported, where as only 779.4 units = 17.4% of the 4475 extra 

wits needed in the second case would have to be imported. 

The over-all cost of the program in the latter case [1, p. 38] would amount to $1097.5 
million, of which $511.3 million is for the processing of high- and medium-grade material into 
alloys. This leaves $586.2 million for the acquisition of that material, which is less than 80% 
of the $744.5 million found in the first computation. The $586.2 million consists of $128.7 mil- 
lion worth of high-grade material contained in the initial stockpile (again evaluated at its 
shadow price), $99.0 million for importation, and $358.5 million for domestic production of the 
additional quantities. From a comparison of these results with those obtained from the first 
computation it is clear that slight improvements in the technology of upgrading and processing 
domestic material would have a great effect not only on the minimum-cost of the program but 
also on the distribution between importation and domestic production of high- and medium- 
grade material. 

Besides these results there were some other interesting features. First, it was 
observed that material flowed into the stockpile before it was completely drained from its 
original contents [1, pp. 18, 27]. This is different from that sometimes found in linear- 
programming problems of this sort. 

Second, and of more interest, was the information contained in the shadow prices inso- 
faras they were related to the constraints. They indicate by how much the total cost of the 
program would rise or fall, should the corresponding restrictions on the variables be more or 
lss severe. They showed that it would be more economical, under the given price and cost 
“nitions, to start the program with higher capacities of some of the plants rather than with a 
larger quantity of high-grade material in the stockpile [1, pp. 24-25, 33-34]. This would indi- 
tate that exercising restraint in the purchase of extra quantities of foreign material would 
te wise, 

Third, there appeared a second set of shadow prices that do not show up in solutions of 
inear-programming problems. These shadow prices indicate by how much the. activities that 
lve zero values in the solution are overpriced. 

All the results described so far were obtained by employing the simplex technique. The 
westion that remained to be answered was: what will the solution be if the negative quadratic 
tms of expressions (22), (23), and (24) are taken into account as well? That question had 
it been answered, despite the two computations that were made, because the unit cost 
ed in the first and those used in the second computation for domestically produced 
uterials were just two rather arbitrarily chosen points of the functions that underlie the cost 
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expressions (22), (23), and (24). The two computations that were made could therefore give 
only a rough idea of the effect that increasing returns to scale in the production of domestic 
material would have on the final outcome. 

However, a new method became available last year that, can handle problems with 
objective functions containing negative as well as positive quadratic terms besides linear 
ones [3]. The computational technique employed in this new method is not the familiar sim- 
plex technique. But, in common with that technique, it first minimizes the objective function 
without its negative quadratic terms. Once the solution of that simpler problem has been 
obtained, then the negative quadratic terms are brought in, after which the search for a lower 
cost-minimum is continued. 

This new method has been used in a third computation. The objective function without 
the negative quadratic terms that was first minimized in this last computation was exactly the 
same as the one minimized in the second computation. In other words, the linear cost of 
producing domestic material in this third computation was also 5 to 10 percent below those of 
the first computation. Once this objective function was minimized, then the negative quadratic 
terms were added to it and the minimizing process continued. The outcome, in units of 1000 
NT of manganese in the form of alloys, was: 


(a) requirements in all 6 years 5330 
(b) contained in initial stockpile 855 
4475 


(c) to be obtained by importation 540.7 
(d) to be produced domestically 3934.3 


In this solution the over-all cost of the program would amount to $1056.0 million, of which 
$495.6 million is for the processing of high- and medium-grade material into alloys. This 
leaves $560.4 million for the acquisition of that material which is $25.8 million below the 
$586.2 million found in the second computation. The value of high-grade material contained in 
the initial stockpile has been set at $128.7 million, the same amount as before. The remaining 
$560.4 - 128.7 = $431.7 million consists of $62.1 million for importation and $369.6 million for 
domestic production of high- and medium-grade material that is required in addition to the 
initial stockpile. 

A comparison of these results with those of the second computation shows a further 
decline in the share that importation has in the acquisition of additional quantities of high-grade 
material from 779.4 to 540.7 units or from 17.4 percent to 12.1 percent of the 4475 units. The 
over-all cost of the program decreased by $1097.5 - 1056.0 = $41.5 million, distributed over 
$511.3 - 495.6 = $15.7 million lower cost of processing, $99.0 - 62.1 = $36.9 million lower 
cost of importation, and $369.6 - $358.5 = $11.1 million higher cost for upgrading domestic 
material. The latter represents an increase of 3 percent of the former cost, while the cor- 
responding quantity increased by 3934.3 - 3695.6 = 238.7 units or approximately 6 percent. All 
this is in agreement with what one would normally expect from taking the negative quadratic 
terms in the cost of domestically produced ores into account. Their effect on the over-all cost 
of the program as well as on the distribution between importation and domestic production 
turned out to be quite significant. 

A final remark might be made as to the character of the various minima obtained by 
these three computations. The first computation produced an absolute minimum of the 
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truncated cost function, i.e., the cost function without its negative quadratic terms; it contained 
the higher linear cost for the upgrading and processing of domestic ores. Given this function, 
the solution is optimal in the sense that there is no other combination of quantities that meets 
the requirements at a lower cost while staying within the given restrictions. The same holds 
for the solution of the second computation; that, too, produced an absolute minimum of a trun- 
cated cost function with somewhat lower linear cost of domestic ores. 

The third computation produced, on the other hand, a local minimum of the complete 
cost function with all its negative quadratic terms included in it. That is, it produced a mini- 
mum in that part of the complete cost function in which the truncated cost function had its 
absolute minimum. But that local minimum will probably not be the absolute minimum of the 
complete function; the latter will in all likelihood be lower than the former. To the extent that 
this absolute minimum has not been found, the problem has not been solved. The method by 
which the absolute minimum of a cost function that contains negative quadratic terms can be 
calculated in a straightforward manner has, however, still to be developed. 
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LINEARITY OF UNRESTRICTEDLY TRANSFERABLE UTILITIES* 
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It is proved that for games with at least three players the assumption 
of unrestrictedly transferable utilities is equivalent to the assumption 
that the individual players' utilities are linear in money. This is not 
the case for two-person games. 











We quote from Ref. [1]: 

"Most of the past work in n-person (game) theory has supposed that, in 

addition to receiving the payoffs prescribed by the rules of the game, the players 

are permitted to make additional transfers—side payments in the delicate lan- 

guage of the theory, considerations or bribes in more direct vocabularies. 

Indeed, a far stronger assumption is made which is generally subsumed under 

the phrase that utility is 'unrestrictedly transferable.' Of course, it is never 

utility as such that is transferred, for utility is a derivative concept, but com- 

modities to which utility can indirectly be attached by the players. To make any 

sense of the elliptic concept of unrestricted transferability and of the mathe- 

matics employed, one must suppose that there exists an infinitely divisible, real 

and desirable commodity (which for all the world behaves like money) such that 

any reapportionment of it among the players results in increments and decre- 

ments of individual utilities which sum to zero according to some specific set of 

utility scales for the players. This can happen if money exists, provided that 

each player's utility for money is linear and that the zero and unit of each utility 

function is so chosen that the conservation of money implies the conservation of 

utility. When else it can realistically happen is obscure." 

It is the purpose of this note to answer the question raised in the last sentence of this 
qotation. We will show that when n = 3 unrestrictedly transferable utilities imply utilities 
breach player that are linear in money. When n = 2 this implication need not hold. 

Let us denote by P the set of possible outcomes p of the game in question, before side 
wments are made. The utility function of each player i is a function of the outcome p and of 
it amount of money x that he gets as a result of side payments after the play is over. We 
ow that this utility function is uniquely determined "up to the choice of a zero and a unit," 
t,up to an additive and a positive multiplicative constant. The assumption that utility is 
westrictedly transferable," as defined in the above quotation, may be formulated as follows: 
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For each player i (i= 1,...,n), it is possible to choose additive and multiplicative constants 
in such a way so that the resulting utility function u; (x, p) satisfies the following conditions for 
each fixed pe P: 


(1) u, is monotonic! in x; 


n = n = 5 : 
(2) Dnt x, = 0 implies rae u; (x; , Pp) nl uj (0,p). 


Mathematically, we are given n real-valued functions u,&,p), ‘etnias u,&,p), where x ranges 
over the real numbers and p ranges over a set P; it is assumed that for each pe P, the u, 


obey (1) and (2). We wish to prove that for n= 3 there are functions c(p) and k; (p) (i=1,...,2) 
defined on P, such that 


(3) u; (x, p) = c(p)x + k; (p) ; 


here c(p) is independent of x and i, and k, (p) is independent of x. 
Fix p. For each i, let £; (x) = fi (x, p) = u; (x, p) - u, (0, p); then f;(0) = 0 and 


n alas n 
(4) Lint % = 0 implies Dink f;(x;) =0. 


Let i and j be distinct players and x an arbitrary real number. If we set x; =X, x; = -x, 
and x, * 0 for m «i, m #j and apply (4), we obtain 


(5) f,(x) = -; (-x) . 


If k =i, k = j, we prove similarly that f,, (x) = -f; (-x). Hence f(x) = f,, (x). Since this holds 
for arbitrary i and k, it follows that there is an f such that 


(6) {(0) = 0 


and f , (x) = fo (x) =...2 f,) = f(x) for all x. From (4) we then obtain that 


n reawre n 
(7) Zj-1 %; = 0 implies Z;_, £(x;) = 0 


and from (5) that 
(8) f(x) = -f(-x) . 
Now for arbitrary x and y, let Xj =X, X9=Y, Xg=-xX-y, and x, = 0 for m >3. From (") 


and (8) it follows that f(x+y) = -f(-x-y) = f(x) + f(y). This result, together with the assumptio 
that the Uj (and therefore also f) are monotonic, implies the linearity of f according to 4 





IThat the utility functions are monotonic is implied by the use of the word "desirable" © 


describe money in the cited passage. The monotonicity assumption may be replaced by th 
apparently weaker assumptionthat there is at least one player whose utility function is bounded 
in some finite interval. Only the most wildly pathological functions fail to satisfy this conditim™ 
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LINEARITY OF UNRESTRICTEDLY TRANSFERABLE UTILITIES 


known result. Hence f;, which is the same as f, is linear in x; in other words, 
u,(s,p) ~ ¥4(0, p) = fj i6s,P) =c(p)x. Setting k;(p) = u;(0,p), we obtain (3). 

We "remark that to obtain (3) it is not necessary to aeoune that (2) holds for all mem- 
bers (Kj, +++» % n) of euclidean n-space. Let us assume that a u. i; » P) = oe u; (0,p) 
in an open sa connected subset V of the hyperplane ZS. =j %;= 0. Then we can deduce (3) for 
all x in the projection 1;(V) of V on the x,-axis.¢ x 

The proof is a slight  emeng of the above proof. Fix p. Let tx, er x) é€V, and 
set f, {@) = u. (e+; »p) - u, i; »p). The proof now goes through as before, with x restricted to 
a small neighborhood of 0; we obtain (3), with x restricted to a small open interval around 
¢. But all of 1 (V ) may be covered by such intervals, and in the places where the intervals 
overlap, we must get the same coefficients c(p) and k, () . A-simple topological argument now 
leads to the desired result. 

If we are willing to add the assumption that the u; are continuous in x, then we can 
extend our result to the case in which V is the closure of an open connected set or lies between 
such a set and its closure. 

The importance of these remarks lies in the fact that, in most practical cases, not all 
possible combinations of side payments actually come into question. For example, a player i 
would not be willing to accept a total side payment x which would yield him a utility ui (x »P) 
smaller than what he can guarantee himself by his own efforts. Thus each player has a certain 
minimum b; below which he will not accept side payments. In the case in which side payments 
are limited by this factor and by this factor only, V is the simplex formed by the intersection 
of the "corner" x, = by perey XZ b, with the hyperplane wa x; = 0. Although V itself is 
not open, in this particular case the linearity of the u; on V follows from their linearity on 
the interior of V, and the continuity of the Uj is not needed. 

It is of interest to ask under what conditions the marginal utilities c(p) are independent 
of the outcomes p. One condition that ensures this is 


(9) u;(0,p,;) = u;&,po) implies u,(y,p,) = u;&+y, Po). 


Intuitively, (9) says that the difference between outcomes Py and Py can be replaced in a 
natural way by the monetary value x. Another case in which we know that c(p) is independent 


of p is that in which P is the cartesian product of sets P, ee 
depends only on the i-th coordinate of p. 

To construct a counter-example for n = 2, it is sufficient to construct Uy and Uy SO that 
fy and fy are nonlinear and f 1) = -f 9(- x). A simple example can be obtained by setting 
Uy (%p) = Uy (x,p) = =x? +k, (), k, ee, arbitrary. To obtain a more “natural” example, we dis- 
tinguish between “pure” meen and “mixed” outcomes, i.e., probability combinations of pure 
outcomes. For pure outcomes p, let k, (P) and ko (P) take integral values only, and set 


n? and the utility u; (x, p) 


seiieatisiisiseeneeeeeeeeecs, 
This holds true even if V depends on 
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(10) uj(x,p) = 27x +ki(p)) + sin 2mx, i=1,2, 


The utilities for mixed outcomes can easily be calculated from (10). This example is more 
“natural” because it satisfies (9). 
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BLEMS 
= It has been suggested that this journal might serve as a medium for publishing problems 
in the area of Logistics — with the idea that other persons might be interested in those prob- 
lems and might submit comments. Readers are invited to submit brief statements on applied 
and theoretical problems in Logistics. Address letters to Managing Editor, Naval Research 
logistics Quarterly, Office of Naval Research, Washington 25, D. C. 


COMMENTS ON 
"A GENERAL THEORY OF MEASUREMENT APPLICATIONS TO UTILITY*! 


The article by Professor Pfanzagl ! represents a significant contribution to measure- 
ment theory in general and to the theory of the measurement of utility in particular. I would 
like to comment briefly on the section of the paper which deals with von Neumann-Morgenstern 


utility, and more explicitly, on the axiom of "consistency" proposed by the author. 
The axiom is given as: 


M(a + b) P(b +c) = M(aPb) +c. 


linterpret its meaning as follows. Suppose that we consider a lottery ticket of the usual form, 
#, for example, ($10, >? $1, >? . The von Neumann-Morgenstern method asserts that, by 
fering various sums of money we can determine some sum of money, let us say $5.50, so that 
the subject is indifferent between this amount and the ticket, that is, so that they have equal 
ttilities. The axiom of consistency appears to assert the following: 


Q) If utility of ($10, + ; $1, 4, and $5.50 are equal, 


then utility of ($12, 7 $3, ay and $7.50 are also equal, 


Were the c from (1) is $2, the proposition to hold for any value of c. 

The axiom seems at first glance to be inherently plausible. It must, for example, hold 
bor anyone who, like the subject in the example above, values alternatives on the basis of com- 
wring expected money values. Even if this is not the case, if, for example, the subject were 
udiferent between the first ticket and $4, it is true that he would pay $6 for the second ticket 


i 
‘Manuscript received June 1, 1960. 
‘hann Pfanzagl, Nav. Res. Log. Quart. 6:283-294 (1959). 
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since he is sure of at least $2 more on it. Professor Pfanzagl shows that if the consistency 
axiom is accepted, then the utility of money function of the subject is determined to within one 
parameter which is fixed by establishing a single point of indifference. The importance of this 
is obvious. The main disadvantage of the von Neumann-Morgenstern method as a device for 
establishing scales of utility in real situations has been the very large number of questions 
which the subject would have to answer. If the consistency axiom can be accepted, then much 
of this difficulty is avoided. 

Both for the reason of applicability and for the reason of acceptability, the axiom should 
not, I believe, be accepted into our thinking. In constructing a scale to measure the utility of 
money, we must determine how the subject feels about $1 million as compared with $10,000, 
$100, $1, and so on. Suppose that a person is indifferent between the ticket ($10,4; $1, }) 
and $4. Let us, for the moment, agree to accept the consistency axiom. Then the subject will 
also be indifferent between the ticket ($1,000,010, 4, $1,000,001, 4) and an amount $1,000,004. 
The question is, what does this tell us about the shape of the utility function in the neighborhood 
of $1,000,000? I suggest that it may very well tell us nothing at all. The subject, whether he 
chooses as a gift the money or the ticket, is absolutely sure in the latter case of receiving at 
least $1,000,000 more than in the first case and will naturally be willing to pay $1,000,000 
more for the ticket since he is sure to get it back anyway. The only real decision, therefore, 
in the second case, is whether he wants the ticket ($10, a $1, >) or $4 sure, which we 
already know from the first comparison. We have no additional real information about how the 
subject feels about having $1 million as compared with having.$1 or $4 or $10. If we really 
wish to determine the entire utility function from the first comparison plus the consistency 
axiom, we can only assume that the slope of the function is everywhere the same as in the 


neighborhood of the first comparison. Reduced to this form, the axiom seems far less 
plausible. 


Next we may question whether or not the axiom might be expected to hold over con- 
siderable amounts, especially if we admit negative payoffs. Suppose that I offer as a gift toa 
person in moderate circumstances a choice of one of the following: $100,000 or a ticket 
($1,000,000, 4; $0, i) . He may very well prefer the sure $100,000, especially if, for exam- 
ple, he needs $50,000 for his wife's operation. Suppose that we subtract an equal amount c, 
say $200,000, from each of the two alternatives. By the axiom of consistency he should still 
prefer the sure amount, which now has become -$100,000. Yet he may very likely choose the 
ticket ($800,000, 4. -$200,000, 4) since this offers him his only chance of escaping from a 
overwhelming burden of debt. 

Professor Pfanzagl makes the following point. If an individual is willing to pay $2 fora 
ticket ($10, $s $1, 5) , it does not mean that the utility of the ticket is equal to the oa, of 
$2 but that the "status quo" has the same utility as a ticket of the form (status quo + $8, 7; 
status quo - $1, 5) since the cost of the ticket must be subtracted from the value of the prize. 
The two will be equivalent only if the axiom of consistency is valid. This is quite true. But 
there is no reason why the two should be equivalent. If we turn our previous example around, 
our subject might have chosen the ticket ($1,000,000, 4. $0, 4 rather than $100,000 had he 
been in comfortable circumstances and wanted a chance at becoming rich: this does not meal 
that he would give you $100,000 for the ticket. It is true that this distinction is not always 
made by writers on utility and that the failure to make it is an unwitting assumption of the 
consistency axiom, but this does not excuse the failure to make the distinction; it is probable 
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that the "lottery ticket" terminology used in the von Neumann-Morgenstern construction has 
confused this point somewhat. 

If we have interpreted the axiom correctly, and if it is true, as Professor Pfanzagl 
indicates, that the axiom has been tacitly assumed in the well-known studies of utility meas- 
urement that he cites, then it is certainly true, as he states, that much of the recent experi- 
mental work on utility measurement will need critical reappraisal, since it seems questionable 
that we shall want to accept the axiom. 


/sgd/ W. G. MELLON 
Econometric Research Program 
Princeton University 





' tot 
ofp 
> prac 
> auth 
Stat 
; tion: 
stud 
16), 


expos 
for L 
autho 


are gi 
the fo 





RECENT PUBLICATIONS 


INTERINDUSTRY ECONOMICS, by Hollis B. Chenery and Paul G. Clark. John Wiley 
and Sons, Inc., New York, 1959. 345 pp., $7.95. 

The contents of this book are precisely indicated by its jacket description, ''A Compre- 
hensive Survey of Both Theories and Practical Applications in the Field of Interindustry 
| Analysis." Interindustry analysis covers all the techniques which can be utilized in order to 
| malyze the interdependence among different industries and also their relationships in regard 
- to the production of goods and services for final demands and allocation of the primary factors 
of production. I would like to emphasize especially the comprehensiveness in the survey of 
practical applications in the book under review. With their rich experience in this field, the 
authors explain the application of interindustry analysis in many countries, including the United 
States. and Italy. A substantial part of the material is taken from the authors' earlier publica- 
tions, which are well known to students in this field. Some of the materials, however, e.g., the 
studies of capital structure at Harvard (pp. 149 - 53) and the U. S. Emergency Model (pp. 271 - 
16), are published for the first time in this book, although their existence has been well known. 
One of the interesting features in the present book is found in Chapter 11, where the authors 
show how the interindustry analysis combined with the linear programming technique can be 
used for the problems concerning economic development plans. 

At present input-output analysis and linear-programming techniques are the main part 
of interindustry analysis. Thus, the survey of theoretical problems in interindustry analysis 
in this book is in a competitive position with Dorfman, Samuelson, and Solow's Linear Pro- 
gramming and Economic Analysis. In the book under review, however, the emphasis is pri- 
marily on the economic aspects of these techniques and less on mathematical formalism than 
inthe book by Dorfman, Samuelson, and Solow. This reviewer believes that this is exactly 
what is wanted in this field because so many recent papers are written with extremely inade- 
quate knowledge of economic facts. 

In the course of a quick reading the present reviewer found one mistake (or incomplete 
exposition). On p. 31 and p. 52, n. 37, a sufficient condition for the existence of the solutions 


or Leontief system is presented. As this reviewer understands the passages there, the 
authors claim that 





(i) ai, = 0, (ii) } aj =1 for all j, and (iii) pa aij <1 for at least one j 
i i 


we Sufficient for the existence of the solution. A counter-example to disprove this claim is 
the following: 
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0.60 0.60 0.45 
A=/0.40 0.40 0.15 
0 0 0.20 


does not converge, and (I - a)7} does not exist. 


(Reviewed by M. Hatanaky) 


THE POLITICAL ECONOMY OF NATIONAL SECURITY, by James R. Schlesinger, 
Praeger, New York, 1960. 292 pp. 

In the words of the author, this book "attempts to present a comprehensive picture of 
the various economic factors bearing on national security."" There is no doubt that a consid- 
eration of this topic, relatively neglected in the outpouring of current writing on problems of 
defense strategy, has long been called for. A book can be fairly reviewed only in the light of 
its purpose and the audience to which it is addressed. Although these are not specifically set 
out in his introduction, the nonacademic phraseology, as well as the origin of the volume ina 
series of lectures delivered by Professor Schlesinger at the Naval War College, suggest that 
the book is intended as a review of the economic aspects of national defense for the literate, 
but noneconomist, reader. In this purpose Professor Schlesinger succeeds admirably. The 
writing is clear and lively, though marred by an occasional cliche, and the explanations are 
clear. The economist reader will find, if nothing new, nothing with which he would be disposed 
to quarrel, and, as Professor Schlesinger points out, "economists are forced to reiterate some 
simple observations, not because they are particularly penetrating, but because they are true." 

After a brief consideration of the role of the economist in the formulation of policy, the 
book considers the question of GNP and other indicators of a nation's capacity for war. The 
historical evolution and the importance of industrial capacity as a decisive factor in warfare 
and the difficulties involved in using and comparing economic indicators are stressed. The 
author then considers the problem of economic mobilization and installation of controls and the 
budget as a tool for allocation of funds. The latter discussion contains a concise survey of 
many of the usually proposed budgetary reforms: in connection with it the reader would do 
well to read the excellent RAND study on Defense Planning and Organization, by Enthoven and 
Rowen. 

A considerable portion of the book is concerned with the international aspects of 
national defense. Three chapters are devoted to the foreign-trade policy of the United States 
and to the economic and political problems of the underdeveloped countries. The author also 
considers the economic growth of the Soviet Union and places considerable stress on the 
sensible, but sometimes neglected, position that we must expect the Soviet economy to grow 
rapidly in the present stage of its development and that we should not be forced into a race t0 
achieve a comparable rate of growth merely for the sake of growth alone. 

The volume would be very suitable for those who, like the reviewer, are involved in 
teaching courses on defense economics where many of the students have only a sketchy eco- 
nomic background, and who have noted the lack of any really suitable basic (up to date) text. 
Professor Schlesinger has considered only the basic economic problems and, in the interest 0 
brevity, has omitted consideration of other economic aspects of defense. In particular, I have 
the feeling that some of the discussion may leave the general reader with the idea that economit 
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theory can't contribute much to the analysis of decision problems beyond the operationally 
feeble advice that we must "balance marginal gain against marginal cost."" But the discussion 
of what we can actually do to solve some of the problems that the author raises belongs, no 


doubt, in another type of book, one which we hope students of the economic aspects of defense 
(like the author) will give us in the near future. 


(Reviewed by W. G. Mellon) 
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